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Abstract 

We present the results of a phenomenological study of unpolarized nuclear structure functions for 
a wide kinematical region of x and Q^. As a basis of our phenomenology we develop a model which 
takes into account a number of different nuclear effects including nuclear shadowing, Fermi motion 
and binding, nuclear pion excess and off-shell correction to bound nucleon structure functions. 
Within this approach we perform a statistical analysis of available data on the ratio of the nuclear 
structure functions F2 for different nuclei in the range from the deuteron to the lead. We express 
the off-shell effect and the effective scattering amplitude describing nuclear shadowing in terms 
of few parameters which are common to all nuclei and have a clear physical interpretation. The 
parameters are then extracted from statistical analysis of data. As a result, we obtain an excellent 
overall agreement between our calculations and data in the entire kinematical region of x and . 
We discuss a number of applications of our model which include the calculation of the deuteron 
structure functions, nuclear valence and sea quark distributions and nuclear structure functions for 
neutrino charged-current scattering. 
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I. INTRODUCTION 



The lepton Deep Inelastic Scattering (DIS) has been since long time a powerful tool 
to probe the structure of hadrons and nuclei at small and intermediate scales. After the 
discovery of the parton structure of nucleons, DIS remains to be the primary source of 
experimental information on the distribution of quark and gluon fields in the nucleon and 
nuclei and a valuable tool to test predictions of QCD. New data from high-intensity electron 
(Jefferson Laboratory) and neutrino (NuMI at Fermilab and JPARC in Japan) beams will 
allow in future to further extend our knowledge of the nucleon and nuclear structure from 
high-precision experiments. 

The role of nuclei in DIS studies is dual. First, it should be noted that the study of 
nuclei at small space-time scales is interesting by itself and it can provide valuable insights 
into the origin of nuclear force and properties of hadrons in nuclear medium. On the other 
hand the nuclear data often serve as the source of information on hadrons otherwise not 
directly accessible. A typical example is the extraction of the neutron structure function 
which is usually obtained from deuterium and proton data in a wide kinematic region. This 
procedure requires, in turn, a detailed knowledge of nuclear effects in order to control the 
corresponding systematic uncertainties. Another example is the determination of nuclear 
parton distribution functions which are universal high-momentum transfer characteristics 
of complex nuclei. Significant nuclear effects were discovered in charged lepton DIS experi- 
ments in 0, 0, B S la H B 113 11) [H; These observations rule out a simple picture 
of a nucleus as a system of quasi-free nucleons and indicate that the nuclear environment 
plays an important role even at energies and momenta much larger than those involved in 
typical nuclear ground state processes. The study of nuclei is therefore directly related to 
the interpretation of high-energy physics from hadron colliders to fixed target experiments. 
The measurements of nucleus-nucleus and proton-nucleus interactions at RHIC [3] and 
LHC will help to clarify the nuclear modifications of the parton distributions, as well 
as to define the initial conditions towards the studies of new states of matter in heavy ion 
collisions. 

The understanding of nuclear effects is particularly relevant for neutrino physics, where 
the tiny cross-section with matter requires the use of heavy nuclear targets in order to 
collect a significant number of interactions. The presence of an axial-vector component in 
the weak current and the quark flavour selection differentiate neutrinos from charged leptons 
and imply a more complex description of nuclear effects in neutrino scattering. The role of 
nuclear corrections to neutrino structure functions has been recently emphasized ^3] after 
the NuTeV collaboration reported a deviation from the Standard Model prediction for the 
value of the weak mixing angle (sin^ 9w) measured in neutrino DIS One of the original 
motivations of the present work is indeed related to the extraction of the weak mixing angle 
from neutrino DIS data of the NOMAD experiment It must be remarked that nuclear 
effects are important not only in the determination of electroweak parameters, but also for 
the understanding of neutrino masses and mixing. The recent high-intensity NuMI and 
JPARC jl^ neutrino facilities offer the possibility to perform a detailed study of nuclear 
effects in neutrino interactions on a relatively short time scale. The construction of a future 
neutrino factory 2^1 would then allow to reach the ultimate precision of the neutrino probe. 

The main experimental information on nuclear structure functions comes from charged- 
lepton scattering DIS experiments performed at CERN 0, 3 li i, Q 0], SLAG 0, IS, 
FNAL |iiE3 and recently at JLab [IllU. The measurements usually refer to the ratio 7^2 
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of the structure function F2 of two nuclei (usually a complex nucleus to deuterium). Addi- 
tional data from the Drell-Yan reaction of protons off nuclear targets are also available 3 • 



From the studies of data on the ratio 7^2 one can separate a few regions of characteristic 
nuclear effects: depletion of nuclear structure functions at small Bjorken x {x < 0.05) known 
as shadowing region; a small enhancement of nuclear structure functions for 0.1 < x < 0.3 
(antishadowing); depletion with a minimum around x = 0.6 -7- 0.7 followed by a rise at large 
X (known as "EMC effect" after the name of the experiment which discovered it). It is 
interesting to note that a clear dependence has been reported only in the shadowing 
region, while for 0.1 < a; < 0.6 TZ2 is almost independent. However, the data available 
on the dependence of nuclear effects are still scarce. One of the main drawbacks of all 
existing data is the strong correlation between x and resulting from the kinematics of 
fixed (stationary) target experiments. As a result, significant regions of the [x, Q"^) plane 
are still uncovered in DIS experiments. 

Many different theoretical models have been proposed to explain the basic features of 
data (for a detailed summary of the current understan ding of nuclear corrections we refer to 
recent reviews and references cited therein (j^. l25l. l26l l27l|. see also discussion in Sec. II VI of 
this paper). The modelling is important to derive some insights on the underlying physics 
of observed phenomena. However, consistent and quantitative description of nuclear effects 
in DIS in a wide kinematical region of x and and for a wide range of nuclei are clearly 
needed. In this paper we perform a quantitative study of data aiming to develop a model 
of nuclear DIS applicable in the analysis of existing data and in the interpretation of future 
experiments. In order to describe nuclear data over a wide kinematical region we take 
into account many effects including nuclear shadowing, nuclear pion excess, Fermi motion, 
nuclear binding and off-shell corrections to bound nucleon structure functions. It should be 
noted that if some effects, such as Fermi motion and nuclear binding, are well constrained 
by other studies or data, the remaining ones are less known. The main example is the off- 
shell correction which describes the modification of structure functions of bound nucleons 
in nuclear environment. We study this effect phenomenologically by parameterizing the off- 
shell correction to the nucleon structure function in terms of a few parameters which are 
fixed from statistical analysis of nuclear data, together with the corresponding uncertainties. 
It is worth to emphasize that these parameters are universal, i.e. common for all nuclei, 
since they are related to the nucleon structure. In a certain sense the off-shell correction 
can be considered as a new structure function which describes the response of the nucleon 
parton distributions to the variation of the nucleon invariant mass. Even if this structure 
function is not accessible for free proton and neutron, it can be probed in nuclear reactions. 

It should be also emphasized that different nuclear effects in different kinematical regions 
of X are correlated by DIS sum rules. For example, the light-cone momentum sum rule links 
bound nucleon and pion contributions to DIS. We use this requirement in order to constrain 
mesonic contributions to nuclear structure functions. Another example is the baryon number 
sum rule which links shadowing and off-shell corrections. In our approach the off-shell effect 
provides the mechanism of cancellation of a negative nuclear-shadowing contribution to the 
normalization of nuclear valence quark distributions. 

After fixing the parameters of our model, we compute predictions for a number of appli- 
cations. In particular, we discuss nuclear valence and sea quarks at high and compute 
nuclear corrections to neutrino structure functions. These subjects will be treated more 
extensively in future publications. 

The paper is organized as follows. In Sec. Owe briefly summarize the DIS kinematics for 
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electron (muon) and neutrino scattering and introduce notations used in this paper. Sec- 
tion mil provides information on the nucleon structure functions and parton distributions 
necessary for our analysis. Section HVl is devoted to the theoretical framework to treat differ- 
ent nuclear corrections in our studies. In particular, in Sec. II V Al we examine the derivation 
of nuclear structure functions in the approximation of incoherent scattering off bound nu- 
cleons and nuclear pions fSec. lIV K^ . the off-shell effects in the structure functions and 
quark distributions fSec. lIV A 6|) and coherent nuclear effects leading to nuclear shadowing 
and antishadowing fSec. lIVB|) . In Sec.|3we discuss in detail the nuclear input which is 
used in our analysis (Sec. lVAl to IV C|) . the model of off-shell effects fSec. lVD]) and effec- 
tive scattering amplitude (Sec. IV E]) . The analysis of data is described in Sec. lVFI and lVG| 
In Sec. lVIl we present the results obtained from our fits to nuclear data. The and A 
dependence of nuclear effects are discussed in Sec. IVI Dl to IVI El In Sec. lVIII we apply our 
approach to study nuclear parton distributions f Sec. lVIl"^ and neutrino structure functions 
(Sec lVIl"^ . In Appendix |X] we provide the details of the integration in nuclear convolution 
and in Appendix ^ the multiple scattering coefficients are given. 



II. KINEMATICS OF LEPTON INELASTIC SCATTERING 

Consider the scattering of a charged lepton (electron or muon) off a nucleon with the 
four-momentum p = {Ep,p) and mass M. The scattering matrix element to leading order 
in the electromagnetic coupling constant a = e^/An ^ 1/137 is determined by the standard 
one-photon exchange process. In inclusive scattering, the final hadronic state is not detected 
and the differential cross section is fully given by hadronic tensor W^iy, which is the sum of 
hadronic matrix elements of the electromagnetic current J^™ over all final hadronic states 
(see, e.g., H) 

W,,{p,q) = ^ Y.^2n)%p + q-pr,) I J7(0)| n) {n iJ^Wlp) , (l) 

n 

where q is four- momentum transfer to the target. 

We do not consider the polarization effects and assume the averaging over the target 
and beam polarizations. Then only the symmetric part of the hadronic tensor contributes 
to the cross section. Because of the conservation of electromagnetic current, time reversal 
invariance and parity conservation in electromagnetic interaction, the symmetric hadronic 
tensor has only 2 independent Lorentz structures (see, e.g., 

F2 

Wi,^{p,q) = -gf,uFi+p^,Py , (2) 

p-q 

where Fi^2 are Lorentz-invariant structure functions, and, for simplicity, we use the following 
notations: 

9t,u = 
h = 

We use the normalization of states {p\p') = 2Ep{27T)^6{p — p') for both bosons and fermions. 
With this normalization the hadronic tensor and the structure functions Fi 2 are dimension- 
less. 



p-q 



(3a) 
(3b) 
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The structure functions are the functions of two independent invariant variables. In deep 
inelastic regime the Bjorken variable x = Q'^/{2p ■ q) and four-momentum transfer squared 

= —q^ are usually used as the variables the structure functions depend on. 

The polarization averaged differential cross section is determined by the structure func- 
tions Fi 2- In terms of the variables x and the cross section reads 



da;dg2 xQ^ 



(4) 



where y = p ■ q/p ■ k. The variable y is not independent variable and related to x and 
via the equation xy = Q'^/{2p ■ k). The lepton mass term is kept in Eq.(j3)) for the sake 
of completeness. Although it is negligible in electron deep inelastic scattering, it might be 
relevant for muon scattering at small momentum transfer or for r lepton scattering. 

The structure functions Fi^2 can be related to the virtual photon helicity cross sections 
by projecting Eq.(^ onto the states with definite photon polarizations. These states are 
described by the photon polation vectors. In the reference frame, in which the momentum 
transfer is along the ^-axis, q = (go, 0_l, g^), qz = —\q\ the photon polarization vectors are 



e± = (0,l,±z,0)/V2, (5a) 
Co = {qz,0±,qo)/Q. (5b) 

where Q = \/Q^. The polarization vectors e+ and e_ describe two transversely polarized 
states with helicities +1 and —1. The vector cq corresponds to the longitudinally polarized 
(scalar) virtual photons. The polarization vectors are orthogonal to momentum transfer, 
e± • g = Co • g = 0, and obey the orthogonality and the normalization conditions, e± ■ Cl = 
0, ■ e± = -1, eg = 1. 

The helicity structure functions are 

W± = e^*W,,el = Fi, (6a) 
Wo = e^W^M-e^ = l'F2/i2x) - F,, (6b) 

where 7 = |qr|/go = (1 + ix^M^ /Q'^Y^'^ . Instead of Eq.®, it is more convenient to use the 
transverse and the longitudinal structure functions defined as 

Ft = x{W+ + W_) = 2xFi, (7a) 
Fl = 2xWo = 72^2 - 2xFi. (7b) 

Let us briefly consider the scattering of (anti) neutrino. In the Standard Model neutrino 
can either couple to charged bosons or to neutral Z boson. In the former case interaction 
is driven by charged curent (CC) = — with and the charged components 
of the vector and axial- vector current. The interaction with Z boson is described by the 
neutral current (NC) which is the superposition of the isovector weak left current and elec- 
tromagnetic current J° = v^(V^ — — 2 siv? 9w J'^^) , where 9w is the Weinberg weak 
mixing angle. 

The hadronic tensor for CC or NC interaction is given by Eq. (^J with the electromagnetic 
current replaced by the corresponding weak current. The Lorentz decomposition of hadronic 
tensor is different for neutrino case and includes additional terms compared to Eq.Q. For 
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example, for CC neutrino interaction we have^ (see, e.g., [2 



p-q 2p-q /g^ 

q^q\j-,w± , q,,Px + qxPf^ j-.w^ 
p-q 

where we denote e^x{ci,h) = e^xapa^b^. The first two terms with Fi and F2 in Eq.Q are 
similar to those in charged-lepton scattering and appear due to VV and AA interactions in 
Eq.(^. The term with F3 describes parity- violating VA and AV interference. The terms F4 
and F5 are present because the axial current does not conserve. The contributions from the 
structure functions F4 and F5 to the neutrino production cross-section are suppessed by a 
small ratio mf/{ME) (these terms vanish in the NC cross section). It was also shown that 
F4 = and 2x^5 = F2 in the leading order and in the limit of massless quarks (Albright- 
Jarlskog relations [l^). Recently it was argued that the second of these relations survives 
the higher order and the target mass corrections in massless QCD, while the relation for F4 
should be replaced by F4 = F2/{2x) - Fi 

The relations between the helicity structure functions and the structure functions -^1,2,3 
in the neutrino scattering are 

W± = F,±^F,, (9a) 
Wo = 7'i^2/(2x)-Fi. (9b) 

The definition of Ft.l by Eq.((7j) also apply in this case. One observes from Eqs.ijH)) that the 
structure function F3 determines the left-right asymmetry in the transverse helicity structure 
functions. 



III. NUCLEON STRUCTURE FUNCTIONS 

The structure functions remain important observables to probe QCD structure of proton 
and neutron and nuclei. In this section we briefiy review the characteristics of nucleon 
structure functions necessary for our analysis. 



A. QCD perturbative regime 

In the region of large compared to the nucleon scale the structure functions can be 
analyzed in perturbative QCD. A working tool of this analysis is the operator product ex- 
pansion (OPE) js^]. Using the OPE, the contributions from different quark-gluon operators 
to hadronic tensor can be ordered according to their twist. For the DIS structure functions 
this leads to the expansion in inverse powers of Q^: 

F^ix, Q') = Ft^{x, Q') + + O (1/g^) , (10) 



The tensor corresponds to interaction mediated by boson and describes neutrino CC scattering 
while describes antineutrino. It should also be remarked that the neutrino and antineutrino NC 
structure functions are identical, since neutrino and antineutrino in NC scattering couple to the same 
hadronic NC. This is not the case for CC neutrino and antineutrino structure functions. 
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where a labels the type of the structure function (a = T, 2, 3). The first term is the leading 
twist (LT) contribution and Ha are the twist-4 contributions (higher twist, HT). 

The leading twist contribution is directly related to the distributions of quarks and glu- 
ons inside the nucleon, the parton distribution functions (PDFs) via the DIS factorization 
theorem as a convolution with coefficient functions (for more detail see, e.g., Ref.js^ and 
references therein). The coefficient functions depend on the process and the type of the 
structure function but are independent of the target. These functions are computable as 
power series in as- The parton distributions are independent of the process but do depend 
on the target. 

The PDFs have non-perturbative origin and cannot be calculated in perturbative QCD. 
However, the dependence of the PDFs can be handled using QCD perturbation theory, 
and is governed by the well-known DGLAP evolution equations with the kernel given by the 
splitting functions [s^. 

The one-loop (NLO) coefficient and splitting functions have been computed since long 
time [s^l- The two-loop (NNLO) coefficient functions [s^ and the corresponding splitting 
functions [s^ are now also available. In our analysis of nuclear data we use both the 
coefficient functions and the PDFs to NNLO approximation calculated in MS scheme using 
the factorization and the renormalization scales set to Q^. 

The HT components involve interactions between quarks and gluons and lack simple 
probabilistic interpretation. 

It must be noted that the twist expansion was derived in the massless limit. If a finite 
mass for the nucleon target is considered, the new terms arise in Eq. (|l(jp that mix operators 
of different spin, leading to additional power terms of kinematical origin - the so-called 
target mass corrections (TMC). If the parameter x^M^/Q^ is small, the TMC series can be 
absorbed in the leading twist term [s^. Therefore, Eq. ()10|) remains valid with the LT terms 
replaced by 

2x^M^ dz 



F™C(x,g^) = -^F^^^(e,g^) + ^^ ^F^^{z,Q'), (11a) 



X 



Fl^'M') = ^^2^"(e,Q^) + ^/ $F-(.,Q^), (lib) 



.„r^ o. 2x^M^ r^dz 



xFj^^ix, Q') = i-.CFt^{^, Q') + ^ / ^zFt^iz, Q% (11c) 



where 7 = (1 -|- 4x^M^/Q^)^/^ and ^ = 2x/ (1 -|- 7) is the Nachtmann variable jS; 



However, it must be remarked that the derivation of |37| was given in the zeroth order 
in as, assuming that the target quarks are on-shell and neglecting the transverse degrees of 
freedom. Furthermore, Eqs. fjllll suffer the so called threshold problem. Indeed, it follows 
from Eqs. (jllll that the target mass corrected inelastic structure functions F™'^ remain finite 
as X — > 1 even if the LT terms vanish in this limit. Clearly, the region x close to 1 is beyond 
the applicability of Eqs. (fTT|) . However, in the applications to nuclear structure functions at 
large x it is important to meet the threshold condition. One possible way to deal with this 
problem is to expand Eqs.(|ll|) in power series in and keep a finite number of terms. In 
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particular, keeping the LT and the 1 /Q^ term we have 



F™''{x,Q')=F^^{x,Q')+ 
x^M^ 









pTMC 


^x, Q 1 


x^M^ 








xF™"" 


^x, Q 1 





(e / Q') - ^i^2'^(^, Q')) , (12b) 

xFl^{x,Q^) + 



2 ^zFl^iz, Q^) - ^ {xFl^{x^ g^))) . (12c) 



In this approximation the structure functions have a correct threshold behavior and vanish 
in the limit of x — >■ 1, provided that the LT terms and their derivatives vanish in this limit. 

The target mass corrections should also be applied to the HT terms in the higher order 
terms in the twist expansion (fTTl|l . For this reason we do not consider 1/Q^ terms in the 
TMC formula, which are small in the considered kinematical range. We also note, that the 
extrapolation of the target mass corrections to off-shell region 7^ is important in the 
treatment of the nuclear effects and will be discussed in Sec. lIV A6l 



B. Structure function phenomenology 

The twist expansion and PDFs as universal, process-independent characteristics of the 
target are at the basis of extensive QCD phenomenology of high-energy processes. In phe- 
nomenological studies, the PDFs are extracted from QCD global fits. A number of such 



analyses are available j39|, |40|, |41|. In our studies of nuclear data described in Sec. lVFI 
to IVI El we use the results by Alekhin js^ who provides the set of the nucleon PDFs ob- 
tained with the coefficient and splitting functions calculated to the NNLO approximation.^ 
Furthermore, the HT terms and the PDF uncertainties have also been evaluated in js^. 

It should be also remarked that the twist expansion and perturbative QCD apparently 
breaks down at low Q^. Furthermore, the conservation of electromagnetic current requires 
the structure function F2 to vanish as for g^ ^ 0. The data seem to indicate the 
presence of a transition region between perturbative and non-perturbative regimes at 
about 1 GeV^. In our studies of nuclear effects in the structure functions some data points 
at small x are in the low-g^ region. In order to match low-g and high-g regions we 
apply spline interpolations for the structure functions which obeys the current conservation 
requirements. 



^ In our analysis we use PDFs obtained from new fits optimized in the low region and including additional 
data with respect to [s^ . This extraction of PDFs also takes into account the nuclear corrections to D 
data described in the present paper (see Sec. IV G|l . Results from the new fits will be reported elsewhere. 
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IV. NUCLEAR STRUCTURE FUNCTIONS 



In this Section we describe a theoretical framework which will be the basis of phenomeno- 
logical studies of nuclear DIS data discussed in Sections |V| to IVll 

The mechanisms of nuclear DIS appear to be different for small and large Bjorken x as 
viewed from the laboratory system. The physics scale for this separation comes from the 
comparision of a characteristic DIS time, which is also known as loffe length Lj = 1/ (Mx) 
(see, e.g., Ref.jl^), and an average distance between bound nucleons in nuclei which is about 
1.5 Fm. At large x > 0.1 the characteristic DIS time is smaller than average internucleon 
distance. This observation justifies the use of the incoherent approximation for the nuclear 
Compton amplitude in this region. It was realized long ago that the nucleon momentum 
distribution {Fermi motion) is important effect even in the scaling limit and results in 
the enhancement of nuclear structure functions at large Bjorken x j43|. After discovery of 
the EMC effect |l|| the calculation of nuclear DIS in impulse approximation was revisited 
"3 S iM S 50l| and effect of nuclear binding was emphasised which explains a significant 



part of the observed dip in the EMC ratio at x ~ 0.6 (for a review of the EMC effect and 



more references see |25l. l26l l27l|). 

Effects beyond the impulse approximation are important. It should be noted that because 
of binding, the nucleons do not carry all of the light-cone momentum of the nucleus and the 
momentum sum rule is violated in the impulse approximation. A natural way to correct 
this problem is to explicitly consider the pion contribution to the structure functions [s^ 
which balances missing momentum. Several calculations of the pion correction to nuclear 
structure functions have been performed in different approaches and approximations js^. 
Although all calculations predict some enhancement at small x, the concrete predictions are 
model-dependent. In this paper we calculate nuclear pion correction following the approach 
of Ref. |47[ in which the pion contribution was constrained using the equations of motion for 
interacting pion-nucleon system. By using the light-cone momentum balance equation we 
effectively constrain the contribution from all mesonic fields responsible for nuclear binding. 

It should be noted that bound nucleons are off-shell particles and their structure functions 
can be different from those of free nucleons. Off-shell effects in nuclear DIS were discussed 



in a number of papers |45|, |46|, |48|, |50|, |5l|, |53|. It was shown that, because of spin 1/2, 
the off-shell nucleon is characterized by the increased number of structure functions which 
depend on the nucleon virtuality as an additional variable [i^ li^ 5^, ^ • However, in the 
vicinity of the mass shell (which is the relevant case for nuclei) the off-shell nucleon can still 
be described by the same number of structure functions as the on-shell nucleon jisl . 49l 50l | . 
Nevertheless, the off-shell dependence of structure functions remains an important effect 
through which the modification of the internal structure of the bound nucleon in nuclear 
environment can be assessed. It should be also emphasized that the off-shell effect provides 
the specific mechanism of balancing a negative contribution to the nuclear baryon number 
sum rule from nuclear shadowing effect (for details see Sections lVl Al and lVl B|) . In Sec. II V Al 
we discuss the derivation of the nuclear structure functions in the presence of off-shell effects 
with the full consideration of the nucleon spin. The treatment of the off-shell effect in the 
parton distributions is discussed in more detail in Sections IIV A 61 and IV Dl 

In the small-x region the space-time picture of DIS is different. For x <C 0.1 the charac- 
teristic DIS time is large on the nuclear scale, the nuclear DIS becomes "stretched" in time 
and in the longitudinal direction. The process can be viewed as the intermediate boson first 
fluctuates into a quark pair which can form a complex configuration (hadronic or quark- 
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gluon) which then scatters off the target. As an average time of hfe of such fluctuation is 
large compared to average distance between bound nucleons, the photon interaction with 
nuclear targets resembles hadronic properties 5^ G^l • In particular, since hadron scattering 
amplitudes are almost imaginary at high energy, the double scattering correction to the 
DIS cross-section is negative leading to nuclear shadowing effect, similar to that in hadron 



scattering [57|. Nuclear shadowing in DIS was subject to intensive studies j6j] (for a review 
of nuclear shadowing and more references see, e.g., (j^. E^). In the present paper we treat 
nuclear shadowing effect in a semi-phenomenological approach by introducing phenomeno- 
logical amplitude which describes interaction of hadronic component of the intermediate 
boson with the nucleon and consider the propagation of this state in nuclear environment 
using multiple scattering theory. Details are discussed in Sec. lIVBl 

Summarizing we write the nuclear structure functions as the sum of incoherent and 
coherent contributions 

where F^^^, Fa^^ , Fl^^ denote the contributions to the structure function of type a from 
bound protons, neutrons, and nuclear pions, respectively. The last term in Eq. (jl3|) is a 
correction due to nuclear coherent interaction. The exact meaning of all these terms will be 
explained in the following sections. 



A. Incoherent scattering approximation 



The DIS hadronic tensor is given by the imaginary part of the virtual photon Compton 
amplitude in the forward direction. In the incoherent scattering regime (large x) ta king into 
account the nucleon spin the nuclear hadronic tensor can be written as (see also 0,13143) 



W. 



=p,n 



tiPA,q) = / [dp] Tr \w;,{p,q)A^{p-A) 



(14) 



where the sum is taken over the bound protons and neutrons, Tr is taken in the nucleon Dirac 
space and the integration is performed over the nucleon four-momentum, [dp] = d'^p/{27i)^. 
In Eq. (ll4p A'^{p]A) is the imaginary part of the proton (r = p) or the neuteron (r = n) 
propagator in the nucleus 



Alf,ip;A) = J dMW^»*-P-(A|x[/J(t,r)vl;^(0)|A) 



(15) 



with "^ait, r) the nucleon field operator and a and (3 the Dirac spinor indeces. The off-shell 
nucleon electromagnetic tensor 'W^p{p, q) is the matrix in the Dirac space. On the mass shell 



p^ -- 

m 



M^, averaging W^i/(p, g) over the nucleon polarizations we obtain the nucleon tensor 



W;Ap,q) = -Ti (^- 



M)wup,q) 



(16) 



In off-shell region, the Lorentz tensor structure of is more involved than the corre- 
sponding structure of the on-shell nucleon tensor. In order to establish the tensor structure of 
we expand the latter in terms of a complete set of Dirac matrices {/, 7", a"^, 7^75, 75}. 
The various coefficients in this expansion must be constructed from the vectors p and g, and 
from the symmetric tensor g^p and the antisymmetric tensor e^^ap. For the symmetric part 
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of we keep only those terms which are even under time-reversal and parity transforma- 
tions, since only such terms can contribute to Fi 2- Keeping only current- conserving terms 
we have 7 independent Lorentz-Dirac structures which can be written as 4^ ^ 




P{t.9u}al'^, (17) 



where g^u and are given by Eq. (jSj) . The curly braces in the last term denote symmetryza- 



tion over Lorentz indeces, i.e. as^hyx = h{anbiy + a^h^). The coefficients f'^''^ in Eq. lfTTj) are 



the dimensionless Lorentz-invariant functions of x, and the nucleon offshellness p^. 

Similar analysis can also be applied to the antisymmetric part of yV^u for the neutrino 
scattering. This term is described by the structure functions -F3 in Eq.®. For off-shell 
nucleon the result can be written as 



^ 2p ■ q 



(18) 





= fi'' 


+ /!" 






F2 


AO) 
— J2 


+ /f 


+ fP - 




F3 


_ AO) 
— /s 


+ /f 


+ /f ' - 





where the coefficients /g''^ are dimensionless Lorentz-invariant functions of x, and p^. 

By substituting Eqs. (fTTj) and (fTH|) into Eq. lfT^ one observes that at = Eqs. ^ 
and (jHI) are recovered with the nucleon structure functions given by 

(19a) 
(19b) 
(19c) 

It should be noted that the Dirac equation is the underlying reason of simplification of the 
Lorentz structure of the hadronic tensor of the on-shell nucleon. 

One important observation which follows from this analysis is that Eq. p4|) does not fac- 
torize into completely separate nuclear and nucleon parts. The off-shell nucleon is described 
by 7 independent structure functions in the symmetric P-even hadronic tensor {f[^^ and 
/g*"*) and 4 independent structure functions f^^^ in the P-odd antisymmetric hadronic ten- 
sor. These functions depend on p"^ as an additional variable and weighted in Eq. (jl4j) with 
generally different nuclear distributions. 

Clearly, the fact that we have to deal with unknown functions not present for the on-shell 
nucleon introduces additional uncertainty in the calculation of nuclear structure functions. 
However, in practice it may be quite sufficient to treat nuclei as nonrelativistic systems. In 
this limit the nuclear hadronic tensor considerably simplifies, as will be discussed in the next 
Section. 



1 . The limit of weak nuclear binding 

Let us now discuss Eq. ()14|) in the limit of weak nuclear binding. We assume that the 
nucleus is a non-relativistic system with small characteristic momentum and energy of bound 
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nucleons, |p| <^ M, \po — M\ M. The antinucleon degrees of freedom are neglected 
in this approximation. A nonrelativistic approximation to Eq. p4|) is derived using the 
relation between the relativistic four-component nucleon field \1/ and the nonrelativistic two- 
component operator ip (for simplicity we suppress the isospin index r) 

where the nucleon operators are taken in a mixed [p, t) representation. The renormalization 
operator 1 — p^/(8M^) is introduced to provide a correct normalization of nonrelativistic 
nucleon field tp, i.e. the operator ip'^i/j is normalized to the nucleon number to order p^/M"^. 

In order to make the nonrelativistic reduction of Eq. p4j) . we separate the nucleon mass 
from the energy pq and write the four- momentum of the bound nucleon as p = (M + e,p). 
We then substitute Eq. fj^UI) into Eq. p4|) and reduce the four-dimensional Dirac basis to 
the two-dimensional spin matrices. In this way we examine all Lorentz-Dirac structures 
in Eqs. ()17|) and ()18p and keep the terms to order e/M and p^ /M"^. The result can be 
summarized as follows: 

^Tr {A{p;A)W,.{p,q)) = -^^^V{e,p)W,.{p,q), (21) 

where Ma is the mass of a nucleus A and 

V{e,p) = Jdt exp{-zet){A\ij\p,t)i:{p,0)\A)/{A\A) (22) 

is the nonrelativistic nuclear spectral function normalized to the number of nucleons in the 
corresponding isospin state 

l[dp]VP'-{e,p) = {Z,N). (23) 

Note that the factor in the left side of (1211) is absorbed in the normalization of nuclear 
states in Eq.(j22I). The hadronic tensor W^^{p,q) in Eq.(pT|) is given by Eq.Q with 

the structure functions 

^2 _ ^(0) A , p2-M2\ (1) p2 



F.{x, Q\p^) = /r ^1 + ) + /r^^ + n", (24a) 

F.ix,Q\/) = /f (1 + ^^) + + /f , (24b) 

F3(a:,g^/) = /f (1 + ^^) + + (24c) 

From Eq. ()2H) we obtain a nonrelativistic approximation to the nuclear hadronic tensor ()14|l 

T=p,n 

which is a basic equation for further analysis of nuclear DIS. 

A few comments are in order. It should be emphasized that the nonrelativistic limit is 
taken with respect to the nucleon momentum. In the derivation of Eq.()21|) we keep terms 
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to order /M"^ and e/M and neglect the higher-order terms. Furthermore, Eq.()2ip is vahd 
for arbitrary momentum transfer q. Note the factorization of the high-energy amphtude 
W^v from the spectral function V which describes the low-energy part of the problem. In 
the vicinity of the mass shell the hadronic tensor involves the same number of independent 
structure functions as on the mass shell. Therefore the problem of "splitting" of structure 
functions in the off-shell region (i.e. the problem of additional nucleon structure functions) 
can be avoided in this region. Equations (|24p give the nucleon structure functions in the 
off-shell region in the vicinity of the mass-shell and it is easy to see that Eqs. ()24|) reduce to 
Eqs. fll9|) at = thus assuring a correct on-shell limit. 

Let us extract the relations between the nuclear and the nucleon structure functions from 
Eq.(j2Hl). Nuclear structure functions are given by Eq.Q with p replaced by Pa and x by 
= Q^/(2iVC4go)- However, it is convenient to consider the nuclear structure functions 
as the functions of the variable x = (5^/(2Mgo) instead of the "natural" nuclear scaling 
variable xa- We then define F^^^{x,Q'^) = F^^^{xa,Q'^) and xF^{x,Q'^) = xaF^{xa,Q'^)- 
In order to separate the structure functions we contract the both sides of Eq. (j25|) with the 
virtual photon polarization vectors © and consider the helicity structure functions. As a 
result we have 

F^ix^Q^) =J2 I [dp]V^ie,p) (l + ^) (f^ + ^-^fA , (26a) 

T=p,n \ ^ / 

Fi{x,Q') =J2f mV^ie^P) (l + ^) (f[ + ^^^j) ^ (26b) 

T=p,n \ ^ / 

where in the integrand F^ with a = T,L,2 are the structure functions of bound proton 
(r = p) and neutron (r = n) with the four-momentum p = {M + e,p), x' = Q'^/{2p-q) = 
x/[l + (5 + 'ypz)/M] is the Bjorken variable for the bound nucleon and p± is the transverse 
component of the nucleon momentum with respect to the momentum transfer. The off-shell 
nucleon structure functions depend on x', momentum transfer square and the virtuality 
= (M-f-e)^— as additional variable. In Eq. ()26p the off-shell transverse and longitudinal 
structure functions are given by equations similar to (|7j) with replaced by p^, i.e. Ft = 
2x'Fi, Fl = 7'^F2 — Ft with 7'^ = 1 + 4a;'^p^/(5^. Using Eqs. pB|) we have for the nuclear 
structure function 

7^F,-(x,g^) = J2 [ [dp]V^is,p) (1 + ^) (7'^ + ^) F;. (27) 

The nuclear structure function F3 can be extracted from the left-right asymmetry in the 
helicity amplitudes, Eq.©. We have jl^ 

xF,^{x,Q') = J2 j\MV\e,p) i^ + ^)^'Fl. (28) 

Equations (j^HI) to allow us to compute the structure functions of a generic nucleus as 
a convolution of nuclear spectral function, which describes the distribution of the bound nu- 
cleons over momentum and separation energy, with the bound proton and neutron structure 
functions. 

We also comment that the transverse motion of the bound nucleon in the target rest frame 
causes the mixture of different structure functions in Eqs. ()26a|) and ()26b|) to order (note 
p\ terms in these equations). This effect on the ratio F^/ F^ was recently discussed in js^ . 
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2. Convolution representation 



If is high enough to neglect power terms in (j^UH^Hj) then these equations can be written 
as two-dimensional convolution. For example, for the structure function F2 we have 

F,^{x, Q') = I dydv [f,,A{y. v)Fi{x/y, Q\ v) + /./^(y, v)F^{x/y, Q\ v)\ , (29) 

Jy>x 

where fp/A{y^ v) and fn/A{y, v) are the proton and the neutron distributions over the fraction 
of light-cone momentum y and the virtuality v = p^. The proton (neutron) distribution 
function is given in terms of the proton (neutron) spectral function as follows (47l . l48l | 

f{y,v) = J[dp]V{e,p)[l + ^)s(^y-l-'-±p-y{v-p'). (30) 

The distribution functions are normalized to the number of bound protons (neutrons), as 
follows from Eq.(j2SI)- Equations similar to hold for other structure functions with the 
same nucleon distribution functions. If we further neglect off-shell effects in the structure 
functions, Eq. (j^^ reduces to the familiar one-dimensional convolution. 

It is instructive to calculate the average nucleon light-cone momentum (y)^ per one 
nucleon. Using Eq. ljHUj) we have 

{y)N = 1 + — — > (31) 

where (e) and (T) are the average nucleon separation and kinetic energies. Because of 
binding effect we have (y)^ < 1 (using our nuclear spectral function from Sec. IV Bl we have 
for iron {y)j^ = 0.966). The missing nuclear light-cone momentum apparently should be 
carried by fields responsible for nuclear binding. In our approach the missing light-cone 
momentum is balanced by nuclear pion field. Note that this situation is qualitatively similar 
to the balance of light-cone momentum in the nucleon in which about a half of the nucleon 
momentum is carried by gluons. However, in nuclear case the fraction of pion light-cone 
momentum is much smaller because the nuclei are weakly-bound systems. The scattering 
from nuclear pions is discussed below. 



3. Pion contribution to nuclear structure functions 

The lepton can scatter off virtual pions which are exchanged by bound nucleons. The 
pion correction to nuclear hadronic tensor can be written as follows {see, e.g., (47| ) 

iy;/^(p^, q) = lj [dk]D^/Aik)w;:,ik, q) (32) 

where W^^{k, q) is hadronic tensor of a pion with four-momentum k and the function -D,r/A 
describes the distribution of pions in a nucleus. The latter can be expressed in terms of the 
pion propagator in a nucleus as 

D^IA^k) = I d'^xexpiik ■ x) {A \(f{x)cp{0)\ A) , (33) 
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where = (y^i, (y92, v^s) is the pion field operator. The vr^ state is described by real pseu- 
doscalar field ^93, while the charged pion states are described by the complex pseudoscalar 
fields {ifi ± Lp2)/V^- The factor 1/2 in Eq.(|H^ is because of the chosen representation of 
the pion field operator in Eq. (j33p in which particle and antiparticle are identical. 

In the further discussion of pion effect it is convenient to consider the normalized pion 
distribution, i.e. independent of normalization of the target state. We define this as follows 

fdtexp(ik„t){A\vj'(k,t)vj(k,0)\A}/{A\A}, (34) 

^(k,t)^j,\eM'k-rMr,t). (35) 

where ip{k,t) is the pion field operator in momentum representation. Using translational 
invariance it is easy to verify that V^^/Aik) = Dj^iA{k) / {2Ma) in the nucleus rest frame. 

In order to extract the pion contribution to nuclear structure functions we contract both 
sides of Eq. ()32j) with the photon polarization vectors. Assuming that the hadronic tensor 
for off-shell pions is given by Eq.(j2)) we obtain from Eq.( 



Q^) = I [dk]V^/Aik) (fco + 7^.) (f^ + ^^^2^) ' (36a) 

Fl^^^ix, Q') = J [dk]V^/A{k) (fco + 7fc.) (fI + ^^F^^ , (36b) 

f F2^/^(x, Q') = J [dk]V^/A{k) {ko + ik,) (y' + i^2^ (36c) 

xF^^'^ix, Q') = j \dk]V^,A{k) {ko + kj^) x'F^, (36d) 

where Fa^^ denotes the pion correction to the nuclear structure function F^. In the inte- 
grand F^ are the structure functions of virtual pion with four-momentum k, kj^ is trans- 
verse component of the pion momentum relative to the direction of momentum transfer and 
x' = Q'^/{2k ■ q) is the pion Bjorken variable. The pion structure functions in Eqs. ()36p 
depend on x', and pion invariant mass k"^ = k"^ — as an additional variable. The 
transverse and longitudinal structure functions are related to Fi and F2 as Ft = 2x'Fi, 
Fl = i F2 - Ft with 7'^ = 1 + Ax'"^ k"^ / Q"^ . The mixture of the structure functions Ft and 
Fl in Eqs. (j36|) is because of transverse motion of nuclear pions, similar to the corresponding 
effect in Eqs. (j26j) for bound nucleons. 

At high Eqs. (l36p can be written in a convolution form. For example, pion correction 
to F2 can be written as 

F^'^ix, Q^) = I dydvU,A{y, v)F^{x/y, Q\ v), (37) 

J x<y 

U/A{y, v) = 2yM j [dk]V^,A{k)5 [y - ^^^^ - k'), (38) 

Similar equations hold for other structure functions in Eqs. ()36p . If one neglects the off- 
shell dependence of the pion structure functions then Eq. (jHTj) reduces to the standard one- 
dimensional convolution with the pion light-cone distribution which is given by Eq. 
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integrated over v. We note that the distribution function by Eq.()38p is antisymmetric 
function, /tt/aI— 1/) = —fw/Aiy)- This property allows us to derive the sum rules for the odd 
moments of the pion distribution function which will be discussed in more detail in Sec. lVCl 



4- Application to the deuteron 

So far the discussion did not refer to particular nuclear target. In this section we apply 
the discussed formalism to the deuteron. The deuteron is an isoscalar bound state of the 
proton and the neutron. The residual nuclear system is, therefore, the proton or neutron 
and the spectral function is given in terms of the deuteron wave function \1'd(p) 

V^'-{e,p) = 21x5 [e-eD + ^ \^ d{p)? , (39) 

where Sd = -2M and p'^/2M are the deuteron binding energy and the spectator nucleon 
recoil energy, respectively. The deuteron structure functions then become 



F^M') = l^,\^n{p)\\l + ^][F^ + ^-^Fi^]. (40a) 

FEM') = I 4^A^D{p)?U + ^-?^](FE + ^-^Fn. (40b) 



(2vr)3 ' V M 

fF,^(x,g^) = /7^|v^D(p)rfl + ^)(7" + ^^)F2", (40c) 



(27r)3 ' V M 



2x'^pI 






2x''pi 






Qx''^pY 


)f^ 





where = (Fp + F^) /2 with a = T, 2, 3 are the structure functions of the isoscalar nucleon. 
The variables of these structure functions are similar to those of Eqs. ()26I27I28|) and we do 
not write them explicitly. 



5. Application to complex nuclei 

Unlike the deuteron, the spectral function of complex nuclei does not reduce to the 
ground state wave function but includes, generally infinite, set of excited residual states 
(this can be seen directly from Eq.(|22j) by inserting the complete set of intermediate states). 
Furthermore, complex nuclei typically have different numbers of protons and neutrons and, 
in contrast to the deuteron case, the calculation of nuclear structure functions requires both 
the isoscalar and the isovector contributions. In order to take into account this effect we 
explicitly separate the isoscalar and the isovector contributions to Eqs. ()26l27l28p . To this 
end we consider generic integrand in the convolution formulas and write 

J2 V^Fl = + pp-«FP-"/2, (41) 

T=p,n 

where we denote = -pp ± -pn ^^^^ pN ^ \_^pp ^ pn^ pp-n ^ pp _ pn ^^le 

structure function of type a. 
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In an isoscalar nucleus with equal number of protons and neutrons Eq. ()4H) is dominated 
by the isoscalar contribution and one generally assumes PP""- = 0. However, it must be 
remarked that this equation is violated by a number of effects even in the isoscalar nucleus. 
The finite difference between the proton and neutron spectral functions is generated by the 
Coulomb interaction and isospin-dependent effects in the nucleon-nucleon interaction. The 
discussion of these effects goes beyond the scope of this paper and we leave them for future 
studies. Instead we focus on the neutron excess effect for heavy nuclei. 

We write the isoscalar and isovector spectral functions in terms of reduced functions Vo 
and Vi as 

pp-n _ {Z-N)Vi. (42b) 

The functions Vq and Vi are normalized to unity as follows from Eq.()23p. These spectral 
functions are quite different. The function Vq involves the averaging over all isoscalar inter- 
mediate states. The function Vi probes the isovector component in a nucleus and its strength 
is peaked about the Fermi surface as argued in Sec. IV Bl The model spectral functios Vq 
and Vi, which are used in this paper, are discussed in Sec iV Bl 

Using Eqs. (j¥T|) and we can write each of the structure function a = T, 2, 3 as 

^/ = ^(Oo + ^(^r">i' (43) 

where the averaging (-Fa)o i denotes the integration in Eqs. ()26|27|28|1 with the reduced spec- 
tral functions Vo and Vi, respectively. 

We conclude this section by commenting that data are sometimes naively corrected for 
the neutron excess effect neglecting Fermi motion and binding effects (as well as any other 
nuclear effects) in the isovector and the isoscalar distributions. As follows from the present 
discussion, the Fermi motion and binding effects are quite different in the isoscalar and the 
isovector distributions in heavy nuclei. If neglected, this effect may cause an additional 
systematic uncertainty in data and a distortion of final results. 



6. Off-shell effects 

The bound proton and neutron are off-mass-shell and their structure functions differ from 
those of the free proton and neutron. The off-shell nucleon structure functions depend on 
the nucleon virtuality p"^ as an additional variable. Therefore, the off-shell effects in the 
structure functions are closely related to the target mass corrections. Target mass effects in 
the off-shell nucleon can be of two different kinds. First, similarly to the on-shell nucleon, 
we have to take into account the kinematical target mass dependence due to the finite p^/Q"^ 
ratio. We assume that this effect is described by Eqs.(fTT|). where the nucleon mass squared 
is replaced by p"^ (this leads in turn to the modification of the parameter 7 and the variable 
^ in the off-shell region). Furthermore, the dependence on p"^ appears already at leading 
twist (LT) at the PDF level as was argued in jiy, 4^, 5^, 5l|. Thus off-shell effects in the LT 



structure functions can be viewed as a measure of the nucleon's modification inside nuclear 
medium. 

Since we treat nuclei as nonrelativistic systems it would be enough to consider the off- 
shell effect as a correction. We expand the nucleon LT structure functions in the vicinity of 
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the mass shell in series in p'^—M^. Keeping only the linear term we have for F2 



Q\p') = F^ix, Q') (1 + 6Mx, Q') ^j^^ , (44) 
xf( n2^ d\nF2{x,Q\p^) 

where the first term is the structure function of the on-mass-shell nucleon and the derivative 
is evaluated at = M^. Similar expressions can be written for the other structure functions. 

The function 5/2 can be related to the corresponding off-shell functions for the nucleon 
parton distributions. The necessary relation can be obtained by writing F2 in terms of a 
convolution of the parton distributions with the corresponding coefficient function according 
to the given order in as- In order to simplify discussion and illustrate this relation we can 
consider the simple leading order expression of F2 



X 



5Ze^(g, + g,), (46) 



where Cj and qi{qi) are the charge and the distribution of (anti) quarks of the type i and the 
sum is taken over different types of quarks. The off-shell function for the parton distribution 
q{x) is defined similarly to Eq. ()45|) . 6fq = dlnq/dlnp'^. Then from Eq. pUl) we have a relation 

F2{x)Sf2ix) =xY,el [q{x)5U{x) + q{x)5f,{x)] . (47) 

One can conclude from Eq. ()47|) that at large x, where the antiquark distributions can be 
neglected, 5/2 is dominated by quarks. For simplicity we neglect the isospin effect and 
assume 5fu = Sfd = Sfg, then 6/2 = Sfg at large x. At small x both, the quark and the 
antiquark contributions, have to be taken into account. 

Off-shell effects in nucleon structure functions were discussed in jil, [s^l using the spectral 
representation of the quark distributions in the nucleon with four-momentum p 

r /'^max 



q{x,p')= ds dtDg/N{s,t,x,p'). (48) 



The integration in Eq. ()48|) is taken over the mass spectrum of spectator states s and the 
quark virtuality t = k"^ with the kinematical maximum tmax = ^^[p^ — s/(l — x)] for the given 
s and p"^. The invariant spectral density -Dg/Tv measures the probability to find in a nucleon 
with momentum p, a quark with light-cone momentum x and virtuality t and the remnant 
system in a state with invariant mass s. 

We conclude from Eq.()48|) that the p^ dependence of quark distributions can have two 
sources: the p"^ term in tmax (kinematical off-shell dependence), and the p"^ dependence of 
the quark spectral function Dg/^ (dynamical off-shell dependence). The kinematical off-shell 
effect causes a negative correction to the bound nucleon structure functions that results in 
an enhanced EMC effect, as first noticed in [i^ 4^. However, if only the kinematical 



off-shell effects are taken into account the number of valence quarks in the nucleon would 
change with p"^. It can be seen directly from Eq.(|48p that the normalization of the quark 
distribution decreases as decreases, provided that the spectral density is positively defined. 
This observation indicates that off-shell effect of dynamical origin must also be present. A 
method to estimate the dynamical off-shell effects minimizing the model dependence was 
suggested in ji^ , in which the conservation of the valence quark number in off-shell nucleon 
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was used as a constraint. A partial cancellation between the kinematical and dynamical 
off-shell effects was found in lis, S^l- However, the off-shell effect in the structure functions 



remains an important correction. In this paper we treat the function 5/2 phenomenologically 
and fix it from nuclear data as discussed in more detail in Section |V1 



B. Coherent nuclear effects 



Nuclear shadowing effect was extensively discussed in the literature. A recent paper 2^ 
provides a review of both data and theoretical models of nuclear shadowing. 

It appears to be a common wisdom that nuclear shadowing is a result of coherent interac- 
tion of hadronic component of virtual photon with target nucleus. The structure functions at 
small X can be presented as a superposition of contributions from different hadronic states. 
We consider the helicity structure functions Wo and W±, as defined in Eq.®, that will allow 
us to discuss nuclear effects in charged-lepton and neutrino interactions on the same ground. 
We have 

WH = J2^,al{s), (49) 

II 

where crj^(s) is the total cross section of scattering of the hadronic state v with the given 
helicity h = 0, ±1 off the target nucleon (or nucleus) with the center-of-mass energy s = 
Q^(l/x— 1) + and the quantities describe the weight of different hadronic states. 

At low the vector meson dominance model (VMP) was proved to be a good tool 
to evaluate nuclear corrections to the structure functions j60|. In this model the structure 
functions are approximated by contributions from a few vector-meson states. The weights 
for the electromagnetic current are Wy = Q"^ / {-w fl){l + Q"^ / m'^)~'^ with the photon-meson 
coupling constants and the vector meson mass. Usually only the lowest mass vector 
mesons (p°, cj, </>) are important at low ^ 1 GeV^. The VMD structure functions have 
strong dependence and decrease as at high Q. In the generalized versions of VMD, 
the higher-mass states including continuum have also been considered that made it possible 
to apply the model at higher Q"^ (see, e.g., [i^]). 

In this paper we approximate the sum over hadronic states in Eq.(|l^ by a factorized 
form 

Wh{x,Q^) = Wh{x,Q'')ah{s). (50) 

where ah is an effective cross section corresponding to helicity h averaged over hadronic 
configurations and Wh is remaining normalization factor. At low the quantity (Jh cor- 
responds to the average over a few vector meson states. As increases, the averaging in 
involves the rising number of active hadronic configurations. Since the relative weight 
of higher-mass states increases with and the cross section decreases with the mass, one 
can qualitatively conclude that an should decrease with Q^. In the approach adopted in this 
paper we will treat ah phenomenologically. 

In this paper we are concerned with the relative effect of nuclear interactions 

^7^,(x, Q\ A/N) = 5W^{x, Q^)/W^{x, Q^), (51) 

where is the nuclear structure function of helicity h subtracted incoherent contribution 
(cf. Eq. (jl3p ). Assuming that the weight factors are not affected by nuclear effects, from 
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Eq.(jHnj we conclude that the relative nuclear correction to the structure functions equals 
the corresponding correction to the effective cross section 

5nH{A/N) = ^cytlal, (52) 

where ba^^ similar to 6W^, is the nuclear cross section subtracted incoherent contribution. 
The problem of calculation of nuclear corrections to structure functions at small x thus 
reduces to the calculation of multiple scattering effects on effective hadronic cross section. 



1. Application to the deuteron 

We now consider this effect in application to the deuteron. In order to calculate the 
shadowing correction we consider hadron elastic scattering amplitude a{s, k) with s the 
center-of-mass energy and k the momentum transfer. We choose the normalization of the 
amplitude such that the optical theorem reads lma(s,0) = a{s)/2 and parametrize the 
scattering amplitude as a = (i + a){a/2) exp^—Bk"^ /2), where the exponent describes the 
dependence on momentum transfer.^ The hadron-deuteron scattering amplitude in forward 
direction can be written as ISj 



= + a" + 5a^ , , , 

(53) 

2 ' 



where and a" are the scattering amplitudes off the proton and the neutron and 5a^ the 
double scattering correction. C2 can be written in terms of the deuteron wave function as 

Cf = ^ 1 d^k^Snik^, fc^)e-^^i, (54) 
Sn{k)= ! d^re'^-^\^D{r)\\ (55) 



Note that Eq. ()53p is the scattering amplidude of an off-shell hadron with four-momentum 
q. For this reason there appears a finite longitudinal momentum transfer ki = Mx{l + 
m^/Q^), which accounts for a finite longitudinal correlation length of a virtual hadron (ki = 
for the scattering of on-shell particles). 

We apply Eqs. and in order to calculate coherent nuclear effects for different 
structure functions. It should be remarked that helicity conserves in multiple scattering 
interactions and the scattering matrix is diagonal in helicity basis. For this reason the 
multiple scattering corrections involve the amplitudes with the same helicity. We also assume 
no isospin effect, i.e. the effective scattering amplitudes of the given helicity are equal for 
the proton and the neutron. Let us first discuss the transverse structure function Ft- The 
relative shadowing correction to the transverse structure function is 



STZrix, Q'\ D/N) = aria^ - 1)^2^/2, (56) 



^ Such dependence is confirmed experimentally and for low mass vector mesons the value of the parameter 
B is between 4 and 10 GeV^ depending on (see, e.g., 0). 
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where ar and are parameters of effective scattering amplitude of transversaly polarized 
virtual photon. A particular model of the scattering amplitude used in our analysis is 
discussed in Section fV El 

Relation similar to Eq. (|56p holds in the longitudinal channel. It follows from Eq. (j5(J|) that 
the ratio R = Fl/Ft equals the corresponding ratio of the cross sections o-L/ax provided 
that the normalization factor w{x, Q"^) is independent of helicity. This holds in the VMD 
for low-mass mesons and we also assume that this is approximately true for the contribution 
from higher-mass states. Thus using Eq. ()56|) and assuming ol = cxt we find that the 
relative shadowing corrections for longitudinal and transverse structure functions is simply 
determined by R 

5n,M\D/N) _ 

57^T(x,Q^D/iV) -^^"^'^ ^' ^^^^ 

where R{x, Q"^) is calculated for the nucleon. 

Eqs. fl56p and (j37j) allow us to compute the nuclear shadowing effect for the structure 
function F2 in terms of the corresponding correction to Ft- Indeed, recalling Eqs.((Z)) we 

1 + -R 

Taking into account (jFfjl we find the factor (1 -|- + R) difference between shadowing 

effect for F2 and F^. 

Let us discuss the shadowing effect for the structure function xF^. This structure function 
is given by the left-right asymmetry in helicity structure functions — W^. Therefore, in 
this case the problem reduces to computing the multiple scattering effect for the difference of 
the corresponding scattering amplitudes. We denote Aa = a^ — a_. The non-zero difference 
Aa is generated because of vector-axial vector current transitions in the hadronic tensor. 
The double scattering correction to Aa can readily be derived from Eq.()53|) 

6Aa^ = 2iAaaTC^, (59) 

where we denote a-r = (a+ + a_)/2. It follows from Eq.(|59p that the relative shadowing effect 
for the cross-section asymmetry is determined by the cross-section (Tt- Using Eqs. (j53p and 
we find 

5nT{D/N) - ^ ^ 

We observe from this equation that the shadowing effect is enhanced for the cross section 
asymmetry by the factor of 2 with respect to the shadoAving effect for the cross section gt 
if we neglect the effect of real part of the amplitudes [63[. To clarify the origin of this 
enhancement we consider a somewhat simplified VMD model with the single vector meson 
(p meson) and the axial- vector meson (ai meson). In this model the structure functions Fj^ 
and Ft in charged-current scattering are determined by the diagonal vector-vector and axial 
vector-axial vector transitions V^A^ — > VN and AN — >■ AN ^ while the structure function F3 
is driven by the off-diagonal transitions V^A^ — > AN and AN VN . The cross section of 
the off-diagonal transitions is much smaller than the cross sections of the direct processes. 
For this reason, XF3 -C F^ at small x. However, the nuclear multiple scattering corrections 
to off-diagonal process V ^ A are determined by a strong cross section of the diagonal 
processes V —>■ V and A A. This becomes clear if we consider the double scattering term 
for the off-diagonal nuclear amplitude. To this order the nuclear scattering proceeds via 
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two steps: the off-diagonal scattering from one nucleon followed by the diagonal scattering 
from the second nucleon. The off-diagonal scattering can be interchanged with the diagonal 
scattering that leads to the factor of 2 enhancement, which appears to have the combinatorial 
origin. 



2. Application to complex nuclei 

We now turn to the discussion of the shadowing effect in complex nuclei. We apply 
the Glauber-Gribov multiple scattering theory to calculate the multiple scattering effect on 
effective cross sections. Let be the nuclear scattering amplitude in forward direction. We 
will assume no isospin dependence of the scattering amplitude, i.e. = a". Then can 
be written as (see, e.g., j60| and references therein) 



= Aa + 5a^, 
(5a^ = m^C^(a), 



(61) 



where a is the corresponding nucleon amplitude and incorporates the multiple scattering 
effects and read as follows 



d^bdzidz2 pA{b, Zi)pA{b, Z2) exp 



21<Z2 



Z2 



dz' {apA{b, z') - ki) 



(62) 



Here pa is the nucleon density distribution normalized to the number of nucleons A and 
the integration is performed along the collision axis, which is chosen to be 2-axis, and over 
the transverse positions of nucleons (impact parameter b). If only the double scattering 
approximation is considered then the exponential factor in Eg . (j62j) should be omitted. The 
exponential factor in Eq.(|H^ accounts for multiple scattering effects (see, e.g., jH^I). 

We note that Eqs. and (jU^ were derived assuming that the wave function factorizes 
into the product of the single particle wave functions and neglecting short-range correlation 
effects between bound nucleons (optical approximation). We comment in this respect that 
the correlations are relevant only if the coherence length Lc = is comparable to the 
short-range repulsive part of the nucleon-nucleon force, which is about 0.5 Fm. This takes 
place at relatively large x, for which shadowing effect is small (see discussion in Ref.|23|). 

The transverse momentum dependence of elastic scattering amplitudes was also neglected, 
since the transverse size of the meson-nucleon amplitude in the impact parameter space is 
of order B~^^'^, much smaller than the radius of the nucleus. 

We first discuss multiple scattering correction to the transverse structure function. The 
relative shadowing correction is determined by effective scattering amplitude of trans- 
versely polarized virtual photon 



5TZt{A/N) = arReii + aTYC^{aT)/2. 



(63) 



If the real part of the amplitude is small then multiple scattering correction is negative be- 
cause of destructive interference of forward scattering amplitudes on the upstream nucleons 
that causes shadowing of virtual hadron interactions. It should be also noted that if the real 
part is large then the interference in the double scattering term is constructive that would 
lead to antishadowing effect. 
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If the coherence length of hadronic fluctuation is small compared to average nuclear 
radius, Lc <C Ra, then the oscillating factor in Eq. ljU^ suppresses multiple scattering effect. 
The onset point of coherent nuclear effects can be estimated by comparing the coherence 
length of hadronic fluctuation L,, with the averaged distance between bound nucleons in the 
nucleus r^N- The coherent nuclear effects take place if the coherence length is large enough 
Lc > f^N'N- Since for any mass of intermediate hadronic state < {Mx)~^ the region of 
coherent nuclear effects is limited to small x for any Q^, x < (MrNN)"^- Nuclear shadowing 
saturates if the coherence length Lc exceeds average nuclear radius that happens at small 
X and the condition Lc ~ Ra defines the transition region with strong x dependence of the 
shadowing correction. 

The rate of multiple scattering interactions is controlled by mean free path of hadronic 
fluctuation in a nucleus (p^cx)"^. If this is small enough compared with nuclear radius, 
which is the case for heavy nuclei, then multiple scattering effects are important. 

It can be easily seen from Eqs. and (j62|) that if the dependence of on the scattering 
amplitude can be neglected, then Eq. (j57j) generalizes to complex nuclei. This corresponds to 
the case when the double scattering saturates the multiple scattering corrections. Generally, 
for heavy nuclei Eq. (j37j) should be replaced by 

6nL{x,Q^A/N) _ Re[{t + aLyC^{aL)] 

STZrix, Q\ A/N) ^^"^^ ^ ^ Re [(z + aryC^iar)] ' ^ ^ 

with ai and qt the effective scattering amplitudes for longitudinally and transversely po- 
larized photons. The relation between the nuclear shadowing effect for F2 and Ft in heavy 
nuclei can be derived from Eqs. lj^ . (jU^ and (jMj) . 

We now discuss the multiple scattering corrections to the right-left asymmetry in the 
helicity scattering amplitudes and the generalization of Eq. ()60|) to heavy nuclei. The multiple 
scattering correction to the difference Aa = follows from Eqs. and ()62|) . can 

be written as 

6Aa^ = % {a\C^{a^) - c?X^{aS)\ , (65) 

where a± are the corresponding nucleon amplitudes. We now use the fact that |Aa| <^ la^l, 
where = |(o+ + O-) is the amplitude averaged over the transverse polarizations of the 
intemediate boson, and expand Eq.(|65|) in Aa keeping only the linear term. We have 

bAa^ = 2iAaaTC^{aT) - AaalC^iar), (66) 



where 



C^(a) = -idC^{a)/da 



d hdzidz2dz^pA{b, zi)pA{b, Z2)pA{b, 2:3) exp 

Zi<Z2<Z3 



i / dz' (a Pa (^7 z') — ki 



(67) 



The first term in the right side of Eq.(|66p is similar to that in Eq . (|59|) . This term is driven 
by the double scattering term (the quadratic a? term in multiple scattering series). However, 
the higher order multiple scattering terms also contribute to (|66p through nonlinear effects in 
C2 and C^. Note in this respect that the expantion of the term in the multiple scattering 
series starts from the tripple scattering term p\. The analytical expressions for and 
calculated for uniform density distribution, which is used in our analysis described in 
Sec. IV Gl are given in Appendix iBl 
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Using Eq. it is straightforward to compute the relative multiple scattering correction 
to the cross section asymmetry 51Zi^{A/ N) = SAa^/Aa that also determines the nuclear 
shadowing effect for the structure function F3. The resulting experession for STZ/^{A/N) is 
somewhat cumbersome in general case and we do not give it explicitly here. It should be 
noted that 6TZ^{A/N) does not depend on the cross section asymmetry Aa for the nucleon 
but does depend on a a = ReAa/ImAa. If we keep only the double scattering term then 
6Tl/^{A/N) is given by Eq. (j6Up also in the case of complex nuclei. However, this relation is 
violated by higher order multiple scattering terms. 



V. DESCRIPTION OF THE MODEL 

In the following we discuss in detail the model which is used to describe nuclear structure 
functions. The model incorporates the treatment of both the coherent and incoherent pro- 
cesses as described in Sec.|TVl We use the model nuclear spectral function calculated in a 
many-body approach f Section IV B|) . The model pion distribution function is constrained by 
light-cone momentum conservation and equations of motion of pion field f Section IV C|) . The 
nuclear shadowing effect is described in terms of effective scattering amplitude of intermedi- 
ate hadronic states of virtual boson off the nucleon (see Sec. llV B|) . In order to describe data 
on nuclear structure functions we explicitly introduce an off-shell correction to the nucleon 
structure functions, which provides a measure of the modification of the nucleon structure 
in the nuclear environment. This effect and the effective scattering amplitude are treated 
phenomenologically in terms of few universal parameters, common for all nuclei, which are 
extracted from nuclear DIS data in a wide kinematic range of x and Q^. The parameteriza- 
tions of the off-shell effect and of the effective amplitude are discussed in Sections IV Dl and 
IV El Sections I V Fl to FV G 31 describe the analysis of data and our main results. 



A. Deuteron wave function 

Nuclear spectral function V describes the probability to find the nucleon with the mo- 
mentum p and the (non-relativistic) energy e in the ground state of the nucleus. We first 
discuss the deuteron for which the spectral function is determined by the wave function (see 
Eq. ()39p ). The deuteron wave function is the superposition of s- and rf-wave states. In the 
momentum space it can be written as follows 



^D,M = yMP) - j Xlm: (68) 

where ipQ and ?/'2 are respectively the s- and d-wave function in the momentum space,^ m 
is the projection of the total angular momentum on the spin quantization axis, Xim is the 
spin 1 wave function with Sz = m, and S^ip) is the tensor operator 

Suip) = 3{ai ■ p){(T2 ■ p) - (Ti ■ (T2, (69) 



^ In terms of the standard wave functions in the coordinate space u(r) and w{r) the functions ■00 and "02 
are ipo{p) = J dr rjo{rp)u{r) and ip2ip) = (S/tt)^/^/ dr rj2{rp)w{r), where jo and j2 are the 

spherical Bessel functions. Note also a different sign of the d-wave term in Ea. (|68|l with respect to the 
wave function in the coordinate space. 



24 



where p = p/\p\ and cri and cr2 are the Pauh matrices acting on the spin variables of the 
bound proton and neutron, respectively. The momentum distribution in the deuteron is 
given by the wave function (jUHj) squared 




V2 4 



) 



(70) 



where the last term in the right side appears due to the tensor operator (|69p . This term 
vanishes after averaging over the deuteron polarizations, which is the case for the present 
paper. The momentum space partial wave functions are normalized according to 



In order to study the sensitivity of our result to the choice of the deuteron wave function 
we used two different deuteron wave functions: the one which corresponds to the Bonn 
potential jH^ and the Paris wave function jH^. These two wave functions have different 
high-momentum component and in this respect represent two extreme situations. 

B. Nuclear spectral function 

The nuclear spectral function can formally be written as a sum over the set of excited 
residual states. This can be seen directly from Eq. ()22|) by inserting the complete set of 
intermediate states. For simplicity we suppress the explicit notations for different isospin 
states and write 



Here the sum runs over the quantum numbers of the states of A—1 nucleons, which include 
the bound states as well as the continuum states, E^~^ and Eq are respectively the energy 
of the residual nucleus (in the recoil nucleus rest frame) and the ground state energy of 
the target nucleus. The residual system balances momentum of the removed nucleon and 
acquires the recoil energy Eji = p^ /2Ma-i- 

The nuclear spectral function determines the rate of nucleon removal reactions such as 
(e, e'p) that makes it possible to extract the spectral function from experimental data. For 
low separation energies (for \e\ < ei;' ~ 30 — 50 MeV) the experimentally observed spectrum 
is similar to that predicted by the mean-field model jo^]. The mean-field model spectral 
function Vmf is given by the wave functions and energies of the occupied levels in the mean 
field jH^. The mean-field picture gives a good approximation to experimentally observed 
spectrum in (e, e'p) reactions in the vicinity of the Fermi level, where the excitation energies 
of the residual nucleus are small jfi^]. As nuclear excitation energy becomes higher the 
mean-field model becomes less accurate. The peaks corresponding to the single-particle 
levels acquire a finite width (fragmentation of deep- hole states). Furthermore, the high- 
energy and high-momentum components of nuclear spectrum can not be described in the 
mean-field model and driven by correlation effects in nuclear ground state as witnessed by 
numerous studies (for a review see jl^). We denote the contribution to the spectral function 
which absorbs the correlation effects as Vcor{£,p)- 




(71) 



V{e,p) = 2nJ2m- l)n, -P mP)\A)f S {e + Ef + Er- E^) . 



(72) 



n 
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In this paper we consider a phenomenological model of the spectral function which incor- 
porates both the single particle nature of the spectrum at low energy and high-energy and 
high-momentum components due to NN-correlations in the ground state. We first discuss 
the isoscalar spectral function which we write as 

Vo{e,p) = VMF{e,p)+V,or{e,p). (73) 

The low-energy part is described by the mean-field spectral function for which we use an 
approximate expression motivated by closure, i.e. the sum over occupied levels is substituted 
by its average value: 

VMF{e,p) = 2nnMF{p)S {e + E^''^ + Er{p)) , (74) 

with Eji{p) the recoil energy of residual nucleus, E^^^ = E^~^ — Eq the nucleon separation 
energy averaged over single-particle levels, and nuF{p) the corresponding part of the nucleon 
momentum distribution. Note that Vmf is not normalized to the number of nucleons since 
a part of the strength is taken by the correlation term Vcov 

The correlated part of the spectral function Vcor is determined by excited states in (f?^ 
with one or more nucleons in the continuum. Following [70] we assume that Vcor at high 
momentum and high separation energy is dominated by ground state configurations with a 
correlated nucleon-nucleon pair and remaining A— 2 nucleons moving with low center-of-mass 
momentum 

\A-l,-p)^i;\p^)\iA-2y,p2)Sip^+p2+p). (75) 

The corresponding matrix element in Eq. ()72|) can thus be parametrized in terms of the 
wave function of the nucleon-nucleon pair embeded into nuclear environment. We assume 
factorization into relative and center-of-mass motion of the pair [t^ 

{{A-2y,p2 mp,)ij{p)\A) ^ C2Ael{k)^CMiPCM)5ipi+P2+p), (76) 

where ipt^ci is the wave function of the relative motion in the nucleon-nucleon pair with 
relative momentum k = (p — Pi)/2 and tpcu is the wave function of center-of-mass (CM) 
motion of the pair in the field of A— 2 nucleons, pcM = Pi + P- The CM wave function 
ipCM generally depends on the quantum numbers of the state of A— 2 nucleons, however the 
corresponding dependence of the ip^-ei is assumed to be weak. Both the wave functions, '0rei 
and ipcM, are assumed to be normalized to unity. The normalization factor C2 describes the 
weight of the two-nucleon correlated part in the full spectral function. 

Using Eq.()76|) we sum over the spectrum of states of A— 2 nucleons and obtain an ap- 
proximate expression for Vcor in terms of convolution of the relative and the CM momentum 
distributions 

Vcor{£, P) =2txCI j dVld^P2nrel(fc)ncM(P2)5(Pl +P2+P) 

5(e+^ + ^^ + E^A . (77) 
V 2M 2Ma-2 ) ^ ' 

Here n^ci and uqm are the relative and the CM momentum distributions, respectively, and 
£■(2) = E^^'^ — Eq is the energy needed to separate two nucleons from the ground state 
averaged over configurations of A— 2 nucleons with low excitation energy. Note that the 
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minimum two-nucleon separation energy 

^(2) = E^-^ - is of order 20 MeV for medium- 
range nuclei like ^^Fe. 

We can further simplify Eq.((77j) if the momentum p is high enough and \p\ ^ |p2|- 
This allows us to take the relative momentum distribution out of the integral over the CM 
momentum at the point k = p. Then we have 

^ 2.Cln,M (. . ^ + . ^ . E'''>) )^^_ . (78) 

where the averaging is done with respect to the CM motion of the correlated pair and M^, = 
M{A — 2)/{A — l) is effective mass of the system of the residual nucleus of A — 2 nucleons and 
the nucleon with momentum pi. In this approximation the high momentum part of nuclear 
momentum distribution is determined by relative momentum distribution in the correlated 
nucleon-nucleon pair embedded into nuclear environment, ncor(p) = C'l'^re^p)- The energy 

spectrum predicted by Eg. ()78p is a broad peak with the maximum at e pV2M and 

effective width \p\Pcm/M with pcu an average CM momentum. 

We perform the averaging over the CM motion of the NN pair in the field of other A— 2 
nucleons assuming that the CM momentum distribution is given by a Gaussian 

ncuiPcM) = (a/vr)^/%xp(-apcM)- (79) 

The parameter a is related to the averaged CM momentum of the nucleon-nucleon pair, a = 
|IPcm)~^- The latter can be estimated from the balance of the overall nucleus momentum 
(70j , {{Yl PiY) = 0, where the sum is taken over all bound nucleons and the expectation value 
is performed with respect to the intrinsic wave function of the nucleus. This gives (Pcm) — 
2{p'^) {A—2) / (A—l) , with (p^) the mean value of the squared single nucleon momentum. We 
consider configurations in which characteristic CM momenta are small. For this reason we 
should also exclude the high-momentum part in estimating (p^) and we will assume that 
this quantity is given by averaging over mean-field configurations. 

Using Eq. (j79p we integrate over the CM momentum in (|78p and the result reads, 

2M 

VcoAe, p) = ncor{p)—VaTT [expi-apf^^J - exp(-ap^^^)] , (80) 

where p = \p\, pmin and Pmax are respectively the minimum and the maximum CM momenta 
allowed by the energy-momentum conservation in Eq. ()77j) for the given e and p, 

Pm.. = M,p/M + pt, (81a) 
Pmin=\M,p/M -pt\, (81b) 

where px = (2M,(|e| - ^th))^^^ and = E'^^^ + Er{p). The latter is the threshold value 
of the nucleon separation energy for discussed configurations. Note that in our notations 
e < 0. We also remark that pt has the meaning of the maximum CM momentum in the 
correlated NN-pair in the direction transverse to p for the given e and p [t^ . 

In numerical evaluations we use the parameterizations for nMF(^, p) and ?7,cor(^, p) of [l^l 
which fit nicely the results of many-body calculation of nuclear momentum distribution. It 
follows from this calculation that low-momentum part incorporates about 80% of the total 
normalization of the spectral function while the other 20% are taken by the high-momentum 
part. The momentum distributions are presented in [t^ for a limited range of nuclei. In 
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order to evaluate the momentum distributions for other values of the nucleus mass number 
A, we interpolate the values of the momentum distributions for each value of momentum 
\p\. For the parameter E^"^^ we take the two-nucleon separation energy, i.e. the difference 
Eq~'^ — Eq between the ground state energies (note that E (2) > 0). The remaining parameter 
of^MF is fixed using the Koltun sum rule [t^ , which is exact relation for nonrelativistic 
systems with two-body forces 

(e) + (T) = 2^5, (82) 

where Sb = E^/A is nuclear binding energy per bound nucleon and (e) and (T) are the 
nucleon separation and kinetic energies averaged with the full spectral function 

{e) = A-' J[dp]V{e,p)e, (83a) 
{T) = A-' l[dp]V{e,p)^. (83b) 

The function Vi describes the isovector component in a nucleus (see Eqs. ()42|) ). In calcu- 
lating Vi we assume that the correlation part of the spectral function Vcor is mainly isoscalar 
and cancels out in the p — n difference. Then V^'^ is determined by the difference of the 
mean-field configurations for protons and neutrons. If we further neglect small differences 
between the energy levels of protons and neutrons then pP"*^ is determined by the difference 
in the occupation numbers of single-particle levels for protons and neutrons. In a complex 
nucleus the deep levels are usually occupied and their contribution cancel out in Vi. The 
Fermi level has a large degeneracy factor and the occupation numbers for protons and neu- 
trons are different. We then conclude that the difference P^""- is saturated by the Fermi 
level and 

V, = \Mp)Me-eF), (84) 

where ep and (pp are the energy and the wave function of the Fermi level (we have neglected 
the recoil of the A—1 nucleus). 

The isovector correction is usually relevant for heavy-weight nuclei in which there exists 
a considerable neutron excess over protons. For such nuclei the Fermi gas model is supposed 
to be a reasonable approximation and we use this model in numerical applications. In this 
model |0f(p)P oc 6{pf — p), where pp is the Fermi momentum which is determined by 
average nucleon density p = 4p|./(67r^). The normalization coefficient can be found from 
condition (|23|) . according to which Vi is normalized to unity. As a result we have 

pfG ^ _ p^^^^^ _ e^)/(4vrp|). (85) 



C. Nuclear pion distribution function 

In calculating the pion effect in nuclear structure functions the relevant quantity is the 
distribution of pion excess in a nucleus since the nucleon pion cloud effect is taken into 
account in the nucleon structure functions. The inspection of Eqs.ipTj) and (PH|) suggests 
that the pion correction is located at small x < pp/M, which is also confirmed by model 
calculations. In this region a good approximation is to neglect jQ^ terms in Eqs. (l36p as 
well as target mass corrections to structure functions. We also assume no off-shell depen- 
dence of pion structure function and base our discussion on convolution approximation by 
Eqs. (P7j) and 
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Before doing model calculations it is important to realize that the pion distribution 
function is constrained by a number of sum rules. The first moment of ()38p gives an average 
pion ligh-cone momentum 

(2/). = / dyyU/A{y) = (^;+> /M, (86) 

where 6'^_^ = [doip)'^ + [dz^p)'^ is the light-cone component of pion energy-momentum tensor. 
The averaging in Eq. ljHUj) means {O) = J d'^r {A\0{r)\ A) / {A\A) for any operator O. It 
is also useful to consider the average which is proportional to averaged over nuclear 
ground state: 

{y-')^ = I dyy-'U/A{y) = M{^') . (87) 

The pion and nucleon fractions of nuclear light-cone momentum are related by the mo- 
mentum balance equation 

{y). + {y)N = ^- (88) 

Equation (jHHj) . although being intuitively obvious, can formally be derived in a meson- 
nucleon field-theoretic model of nuclear Hamiltonian ji^]. Several constraints on nuclear 
pion distribution ■D^/^(A;) can be obtained in this model using the equations of motion for 
pion and nucleon operators and energy-momentum conservation condition. In particular, 
for a model nuclear Hamiltonian with nucleons and pions with pseudo-scalar interaction we 
obtain the following relations [i^] 



m 



2 



y) =eB + (T) , (89a) 
{ido^)')=eB-'^{{e) + {T)), (89b) 
((V^)^> = -|(5)-i(r). (89c) 

A few comments are in order. Pion field in nuclei is mainly generated by nucleon sources. 
Time variation of the pion field describes retardation effects in the nucleon-nucleon interac- 
tion. In a nonrelativistic system this effect is small since typical energy variations are small 
compared to the pion mass. We, therefore, take the static approximation dQ(f = 0. Then 
Eq. ()89b|) is equivalent to the Koltun sum rule (|82j). In the static approximation for the pion 
energy-momentum tensor we have (6*++) = | ((V(/^)^). Then using Eq. (p?T|l we conclude 
that Eqs. ()89cj) and (jHH|l are equivalent. For this reason only Eq. ()89a|l gives independent 
constraint. 

We use the constraints on average pion light-cone momentum y and 1/y which follow from 
Eqs. fl89a|) and (jHHj) in order to evaluate the pion contribution to nuclear structure functions. 
It should be remarked that in this approach by using momentum balance equation (jHHj) we 
effectively take into account the contributions from all mesons. In order to quantitatively 
evaluate the pion effect in the structure functions we use a model distribution 

U/A{y)=Cy{l-yr, (90) 

which is motivated by the asymptotics of pion distribution function at small and large y. 
The normalization constant C and the exponent n are fixed from Eqs. ()88p and (jHTj) using 
Eqs. fl31|) and ()89a|) . The nucleon average separation and kinetic energies are calculated 
with the spectral function described in Sec. lVBI 
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D. Parameterization of off-shell effects 

The off-shell effect in the structure function F2 is described by Eq. fl45|) . In the analysis of 
data, described in detail in Sec. lV Gl we consider a phenomenological model of the off-shell 
function 5/2 (x, Q"^). In order to choose an appropriate model we first note that function (j45p 
describes the relative off-shell effect on the LT structure function and we expect that Q^- 
independent 5/2 is a good approximation. We also note that off-shell effects are constrained 
by the normalization of nuclear valence quark distribution (see Sec. lVI A]) . For this reason we 
anticipate that Sf2{x) should have at least one zero. Moreover, the analysis of nuclear pion 
correction as discussed in Sec. lV CI suggest that 5/2 (x) can have two zeros. These motivate 
us to choose the following simple parameterization for the off-shell function: 

Sf2{x) = Cn{x- xi){x - Xo){h- x), (91) 

where Cat is an overall normalization constant and < xi < xq < 1 and h > 1. The analysis 
of data indicates that the parameters h and Xq are fully correlated and suggests h = 1 + xq. 
After imposing such condition then expression (j91|) has only three independent parameters. 
We use this model to describe off-shell effects in the analysis of Sec. lV G'2\ 



E. Effective scattering amplitude 

As discussed in Sec. lIVBl the coherent multiple-scattering nuclear effects are determined 
by effective (averaged over hadronic configurations of the intermediate boson) scattering 
amplidudes ah for different hehcities h = ±1,0. The amplitudes at = o"/i(^ + o^h)/'^ are 
parametrized in terms of effective cross section a and the Re/Im ratio a. For the unpolarized 
case, which is considered in this paper, the relevant amplitudes are the average transverse 
Qt = (a+ + a-)/2 and the longitudinal ql amplitudes. One can qualitatively expect that 
ax descreases with since the relative weight of higher mass states increases with and 
the cross-section decreases with the meson mass (see Sec lIVB]) . In order to parametrize 
effective transverse cross section we use then the following expression 

= ^ (i + QV«- 

The parameter ctq describes the cross section at small Q^, while ai corresponds to high- 
q2 YQg[jYie, The choice of both these parameters will be discussed in detail in Sec. lV (ill 
The free parameter describes the transition between low- and high-Q^ regions. We note 
that in the discussed approach we only consider relative corrections to the effective cross- 
section and for this reason the analysis is not very sensitive to the detailed modelling of 
such cross-section. The presence of non-zero real part of the amplitude is required by both 
theoretical arguments and phenomenology. The choice of ar in our analysis will be discussed 
in Sec. lVTTTl 

In order to fix the effective amplitude in the longitudinal channel we use the relation 
(It Jar = R = Fl/Ft with R calculated using the PDFs and the structure functions of 
Ref. j39[ as discussed in Sec. lllll fsee also Sec. II V Bl) . 

The C-odd asymmetry Aa = a+ — a_ in the scattering amplitude of left- and right- 
polarized virtual boson does not affect the structure functions Fi and F2. However, Aa 
is relevant for F3 and affects the normalization of nuclear valence number as described in 
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Sec. IVI Al It should be noted that the relative nuclear shadowing correction to F3 does 
not depend on the cross-section asymmetry Aa but does depend on a/\ = ReAa/ImAa 
as explained in Sec. lIVB'2l In order to fix q;a we use the approach based on Regge phe- 
nomenology of high-energy hadronic amplitudes and approximate Aa by the cj-reggeon pole, 
a simple proper contribution to the C-odd amplitude. The energy dependence and Re/Im 
ratio of the Regge pole is fully determined by its intercept which is about 0.5 that leads to 
«A = 1 IZ^l- We use this value in the calculation of nuclear shadowing correction to the 
valence quark distribution in Sec. IVI Al 

F. Nuclear data 

Table H] summarizes the list of experimental data used in this paper. They include both 
muon (EMC, NMC, BCDMS, FNAL E665) and electron (SLAG E139, E140) scattering on 
a variety of targets: p, D, ^He, ^Li, ^Be, ^^C, 27A1, ^oca, ^^Fe, e^Cu, ^^^Ag, ^^^Sn, ^^^Au, 
^'^^Pb. For each target and kinematic region, we select the most precise and recent data and 
we do not use earlier results characterized by larger uncertainties, since their contribution 
to the present analysis would be negligible.^ Most of the data come from NMC for the small 
X region and SLAC E139 for the region x > 0.1. 

We note that, since all available nuclear data are provided by fixed-target experiments, 
there is always an implicit correlation between x and in data points. Usually low- a; 
regions also correspond to low-Q^ values. As described in the following, this reduces the 
possibility to test the dependence of the model in a complete way. 

G. Extraction of parameters 

The numerical values of the parameters in the model are determined from the data listed 
in Table U with two main steps. Initially, we verify the consistency of our model with F2 
data from charged-lepton scattering, without imposing specific constraints. We then discuss 
in detail the deconvolution of different physical effects which contribute to the overall nuclear 
modification of the structure functions. 

It must be noted that the extraction of parameters responsible for nuclear effects is 
correlated with the determination of PDFs (see Sec. lIIl B|) . which requires both the proton 
and deuterium data to obtain the distributions of d and u quarks. Nuclear effects can produce 
significant distortions on parton densities and will be the subject of a future publication. In 
principle, using our approach it would be possible to extract simultaneously the proton PDFs 
and the parameters responsible for nuclear effects (such as off-shell correction and effective 
cross section) by applying QCD analysis to the extended set of data including nuclear data. 
However, in order to limit correlations we follow a different approach. 

The parameters of the model are extracted only from the measured ratios TZ2{A'/A) = 
F2 jF-ti where A and A are two different nuclei (usually the denominator corresponds to 
deuterium). The description of the nucleon structure functions largely cancels in the ratios, 
thus effectively removing the correlation with PDFs. In order to verify this we applied an 



^ Note also that the addition of unnecessary data points with large uncertainties can produce an artificial 
reduction of the of fits. 

^ We note that D/p data were not used in our fits. We compare our predictions with these data in Sec. IVI El 
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TABLE I: List of nuclear data used in the analysis. See text for details. 



iterative procedure. We first extracted the parameters using PDFs obtained without our 
nuclear corrections. Then we repeated the procedure after updating the PDF extraction 
using the information on nuclear effects in deuterium from the previous step f Section IIIIB|) . 
Results indicated that the fitted parameters were stable, demonstrating the absence of strong 
correlations. 

Nuclear data are usually available in bins of x (Ax), while only the average in each bin 
is provided. We perform a fit to the experimental data with MINUIT [t^ by minimizing = 
J2(J^T^ / ^'^{T^T^)^ where a'^(JZ'^2^) represents the uncertainty on the measurements 

and the sum includes all data points. For each experimental point, the model is evaluated 
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at the given average Q and integrated over the size of the x bin: 

Both the normahzation and point-to-point uncorrelated uncertainties, as pubhshed by ex- 
periments, are taken into account. We would hke to emphasize that the lack of knowledge 
of the experimental distribution in the x bins can potentially result in a mismatch be- 
tween data and predictions in the regions where a significant dependence is present. As 
discussed in the following, this increases the systematic uncertainties of the calculation from 
the measured parameters at x > 0.70 and x < 0.05. 

As explained in Sec. IIII Bl we use a phenomenological extrapolation of free nucleon struc- 
ture functions for < 1.0 GeV^ and, in general, nuclear corrections to structure functions 
can be calculated at low Q^. However, we restrict the fits to extract the free parameters 
of our model to the data with > 1.0 GeV^ in order to reduce systematic uncertainties 
on the parameters. We then validate our predictions against the data points with < 1.0 
GeV^, which are included in all comparisons shown in the following. 



1. Choice of fixed parameters 

We start our fits by treating ctq and ai, the asymptotic values of the effective trans- 
verse cross-section in Eq. (j92p . and the real part of the effective scattering amplitude 
ax = ReaT/Imay f Section IIV B 1|) as additional free parameters. This procedure allows a 
preliminary estimate of their correlation with the remaining parameters and a consistency 
check with the expected values. 

The best fit value obtained for ai is consistent with zero. By setting cri 7^ we can 
still obtain an acceptable description of data provided (Ji < 1.0 mb (at 90%CL), due to the 
(anti)correlation of ai with Qq in Eq. ()92|) . After verifying that the correlation between ai 
and the off-shell parameters in Eq. ()91|) is negligible, we then fix ai = in all our fits. 

We note that a non- vanishing shadowing correction at large affects the normalization 
of the valence quark number per nucleon (Section IVl A|) . In this respect we have two 
possible constraints at large Q^. The first condition is to require the conservation of the 
overall valence number in nuclei through a balance between the shadowing and the off- 
shell corrections. As a second choice, it is also possible to explicitly impose the off-shell 
effect to conserve the valence quark number of the off-shell nucleon. This implies that both 
the off-shell and the shadowing effects conserve independently the normalization of valence 
quark distribution. In our approach, initially we do not assume any specific normalization 
constraint. Instead, we verify a posteriori the magnitude of the renormalization introduced 
by the off-shell effect fSection lVI A|) and its balance with the shadowing correction. We then 
use the normalization condition for nuclear valence number to further bound some of the 
parameters. This procedure will be discussed in more detail in Sections IVTT2I and NTM 

The fits to DIS data on nuclear targets show a strong (positive) correlation between (Tq 
and ax- In addition, the value of ao is also correlated with so that it is not possible 



^ In a few cases, in which the expUcit dependence is provided, the model is averaged over the corre- 
sponding bins. 
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to unambiguously disentangle the three parameters from the fits. If we fix = we 
obtain (Tq = 36 mb from data. However, data clearly prefer ar 7^ 0, with a somewhat 
lower value of ctq. The best fit solution corresponds to ar = —0.179 ± 0.038(stat.) and 
Ax^ ~ 29 with respect to the fit with fixed = 0. This can be interpreted as the evidence 
for a sizeable real part in the effective scattering amplitude. If we impose o"o = 27 mb, 
as expected for electromagnetic interactions by averaging lj and vector mesons, we 
obtain = —0.182 ± 0.037(stat.). Note that this is in a good agreement with the analysis 
of photoproduction experiments j6l| at low-Q^. Since we require our phenomenological 
model to correctly reproduce the photoproduction limit, we fix ctq = 27 mb and ax = —0.20 
according to |6ll |. 

In our model we use the pionic parton distributions extracted from real pion scattering 
data ji^] to approximate the structure functions of virtual pions in nuclei. To this end we 
perform fits with and without the pionic sea distributions and we find a significantly better 
description of data in the latter case. Therefore we only consider the valence contribution 
to the pionic structure functions in the following. 



2. Results 

In our model we assume three main free parameters: C^, Xq and (see Sections lV DI and 
IV E|) . In addition, the off-shell function Sf2{x) is characterized by the presence of a second 
zero, Xi. This specific feature has important consequences, as it is discussed in Sec. lVl iJl 
Since the parameter Xi turns out to be strongly correlated with Cn and Ql, we perform 
several fits with different fixed values of xi in the range 0.030 < xi < 0.065 and we evaluate 
the corresponding effect on the normalization of the valence quarks at large Q^. Among 
the fits with comparable with respect to data, we choose a fixed value xi = 0.050, since 
this value provides a good cancellation between off-shell and shadowing corrections in the 
normalization for all nuclei (see Sec. lVl Al for details). 

In order to test our hypothesis about the universality of parameters in Eqs. (|91|) and 
(I92p . we perform independent fits to different sub-sets of nuclei (from ^He to ^°''Pb) and we 
compare the corresponding values of the parameters with the ones obtained from a combined 
fit to all data. As can be seen from Table ITTl the results are compatible within uncertainties, 
thus allowing a unified treatment. The values of x^/d.o.f. indicate an excellent consistency 
between the model and the data points for all nuclei. 

The final values oi Cn, xq and Ql obtained from a global fit to nuclear data are given in 
the last line of Table ITTl The correlation between the parameters is small and mainly related 
to the normalization constant, as can be seen from Table ITTTl 

Figures Q and |21 show the excellent overall agreement between the calculation and the 
data points for many different nuclei. A few comments are in order. The region at x > 
0.75 is characterized by a significant dependence and therefore the calculation based 
upon the average provided by the experiments is approximate. It must also be noted 
that in some cases the data points from different experiments are not fully consistent. In 
particular, the data points on ^^C/D and '*°Ca/D ratios from E665 experiment 3| at low x 



® Unfortunately it is not possible to have data points covering both the high and low x regions for all nuclei. 
This can result in a small sensitivity to some of the parameters for specific nuclei, as can be seen from 
the uncertainties in Table ITU 
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27A1/D; 27A1/12C 


7 39 ± 86 


470 ± 016 


3.60 ± 2.41 


15.1 / 33 


40Ca/D; ^OCa/i^C 


7.09 ± 0.79 


0.482 ± 0.016 


1.67 ±0.15 


56.8 / 58 




8.28 ±0.53 


0.449 ± 0.009 


1.39 ±0.30 


55.6 / 63 


i08Ag/z); ii9Sn/i2c 


9.61 ± 1.29 


0.448 ± 0.018 


1.41 ±0.34 


21.0 / 29 


i97Au/Z);207pb/D; 207pb/i2c 


8.52 ±0.87 


0.387 ±0.026 


1.31 ±0.37 


18.2 / 42 


All data 


8.10 ±0.30 


0.448 ± 0.005 


1.43 ±0.06 


458.9 / 556 



TABLE II: Values of the parameters extracted from independent fits to different sub-sets of data 
with > 1.0 GeV^. Uncertainties are startistical only. The column on the right gives the 
from each fit and the corresponding number of degrees of freedom. The last row shows the result 



Cij 


j = CN 


J = 


3 = Ql 


i = Cn 


1.000 


-0.067 


-0.127 


i = xq 


-0.067 


1.000 


0.006 


i = Ql 


-0.127 


0.006 


1.000 



TABLE III: Correlation coefficients between the nuclear parameters from the global combined fit. 



seem to be systematically above the corresponding NMC measurements, which have smaller 
uncertainties. Similarly, a normalization problem could be present for ^°^Pb/D data from 
E665. Assuming the effect is common to all heavy targets, in our fits we use instead the 
double ratios (^°Ca/D)/(i2C/D) and {^^'^Vh/T>) / {^"^C/T)) and the E665 measurement of the 
ratio ^°^Pb/D. The double ratios are in good agreement with NMC data (noticed also in j^) 
as well as with our predictions, while the ^°^Pb/D points lie slightly above our calculations. 
Futhermore, the ratio ''Li/D shown in Fig. Q indicates a small excess in the region of x 
between 0.01 and 0.03, which produces corresponding reductions in the ratios ^^C/^Li and 
^°Ca/^Li. The effect is much larger than the quoted systematic uncertainties. For instance, 
the exclusion of three points at x = 0.0125, 0.0175, 0.0250 from our fit leads to the reduction 
of overall x^/d.o.f. for the ratio ^Li/D from 1.95 to 0.72.^ We also comment that the value 
of the ^'^''Pb/^^C ratio at x = 0.7 from NMC (Fig. |2)) is marginally compatible with the 
corresponding value of ^^''Au/D ratio from E139 experiment. 

It should be emphasized that at low x there is an interplay between the off-shell function, 
the pion contribution and the coherent nuclear effects. This results in significant correlations 
between the corresponding parameterizations and does not allow an unambiguous extraction 
of individual components without external constraints. In our approach the pion (meson) 
excess in nuclei is calculated as described in Sec. IV CI In order to disentangle the actual off- 
shell function from the remaining coherent correction, we use additional information from 



^ The comment is only intended to quantify the effect. Wc keep all data in our fits, regardless of the 
inconsistencies described above. 
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FIG. 1: Ratios TZ2ix,A'/A) for ^He/D, ^^C/D, ""^Ca/D, "^Li/D and ^Be/D (left to right and top to bottom). The curves with open squares 
show the corresponding model calculations with the parameters specified in the last line of Table ITl| For data points the error bars correspond 
to the sum in quadrature of statistical and systematic uncertainties, while the normalization uncertainty is not shown. 
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FIG. 2: Ratios 'R2{x,A'/A) for ^^Fe/D, ^^Cu/D, i^'^Au/D, ^o^pb/D, ^OCa/^^Q g^j^^j 207p|3/i2Q ^^f^ ^.-gj^^ ^^^^^ ^^p bottom). The curves 
with open squares show the corresponding model calculations with the parameters specified in the last line of Table ITll For data points the 
error bars correspond to the sum in quadrature of statistical and systematic uncertainties, while the normalization uncertainty is not shown. 



photoproduction experiments (Section IV G The agreement between our independent 
extraction of the average VMD parameters and the photoproduction hmit makes us confident 
in the deconvohition of different components. 

We further check the interplay between nuclear pion excess and off-shell effects in our 
analysis by fitting our model without pion correction to 7^.2 data. In this case the effective 
off-shell function 5/2 absorbs the nuclear pion contribution to nuclear F2. We use for this 
test a higher order polynomial with respect to Eq.(|91|). without any fixed parameter. This 
is intended to avoid biases from the functional form used to model the off-shell function. 
The results obtained for the effective 5/2 are consistent with the following estimate which 
can be obtained by explicitly separating the nuclear pion contribution to nuclear F2 



{v) F^^ [x) 



where SF2^^ is the nuclear pion correction calculated as described in Sec. IV CI and (v) de- 
notes the nucleon virtuality f = (p^ — M'^)/M'^ averaged over the nuclear spectral function. 
Moreover, the best fit corresponds to a value of Qq which is in agreement with our fit with 
explicit treatment of nuclear pion correction. 

Figure El illustrates different nuclear corrections to the ratio of F2 of ^^^Au to that of 
the isoscalar nucleon F^ = + F2) calculated in our model using the final parameters 
shown in the last line of Table HTl As follows from comparison of Fig.lHland the results of our 
fit displayed in Figs. C] and El the standard Fermi motion and nuclear binding effect treated 
in impulse approximation does not quantitatively explain data at large x. The off-shell 
effect is therefore an important correction which modifies the structure functions of bound 
nucleon and affects the slope and the magnitude of the ratio 7^.2 at large x. As discussed 
above, we extract this correction from inclusive nuclear DIS data. The disscussion of off-shell 
correction in terms of a scale characterizing valence quark distribution and its modification 
in nuclear environment is presented in Sec. I VI (!1 



3. Systematic uncertainties 

Systematic uncertainties of the model are evaluated by varying each of the contribu- 
tions from the deuterium wave function, the spectral function, the parton distributions, 
the pion structure function and the functional forms of Sf2{x) and ctt(Q^). New fits are 
then performed and systematics are defined from the corresponding variation of the nuclear 
parameters and from the global values. Results are listed in Table IIVI 

Although we do not use directly deuterium data for the fits, most of the data points 
come from the ratios 7I2 of a heavy target to deuterium. In order to study the sensitivity of 
our result to the choice of the deuteron wave function we performed independent fits with 
two different choices of the deuteron wave function: the one which corresponds to the Bonn 



potential j65| and the Paris wave function j66|. These two wave functions have different 
high-momentum component and in this respect represent two extreme situations. 

Similarly, we modify the high-momentum component of the momentum distribution 
ricorip) in nuclei by multiplying it by the ratio of the Bonn and Paris deuteron wave functions 
squared. This is motivated by the observation 7^ 71 1 that the momentum distribution of fi- 



nite nuclei and nuclear matter at high momenta are proportional to that of the deuteron. We 
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FIG. 3: Different nuclear effects on the ratio of ^^''Au to isoscalar nucleon for F2 at = 10 GeV^. 
Tlie labels on the curves correspond to effects due to Fermi motion and nuclear binding (FMB), 
off-shell correction (OS), nuclear pion excess (PI) and coherent nuclear processes (NS). Target mass 
and the neutron excess corrections are included. 



Systematics 


6Cn 


5x0 


5Ql [ GeV^] 


Deuteron wave function 


0.238 


0.0016 


0.026 


Nuclear spectral function 


0.451 


0.0046 


0.021 


Parton distributions 


0.005 


0.0007 


0.023 


Pion structure function 


0.050 


0.0020 


0.065 


Functional form 


0.120 


0.0040 


0.070 


Cross-sections ctq ^-nd ai 


0.015 


0.0005 


0.165 


Total 


0.526 


0.0067 


0.195 



TABLE IV: The estimate of systematic uncertainties on the extraction of nuclear parameters. See 
text for a detailed explanation. 



then repeat our fits with modified spectral functions in order to estimate the corresponding 
variations on the parameters of the model. 

The systematic uncertainty related to the parton distributions is estimated by varying 
the PDFs within their uncertainty (ila). In addition, we also use different sets of parton 
distributions, extracted from fits to different data samples, with different boundaries and 
different parameterization for the low extrapolation. 

For the pion structure function, we repeat our fits by using both the LO and the NLO 
approximations of the pionic parton distributions from j42|. We also arbitrarily change the 
parameterizations 4^1 within ±10%. 

We tried different functional forms in Eqs. (|91|) and (j92p to parameterize 5f2 and (Jt- 
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In particular, for the off-shell function 5f2 we have tried a generic higher order polinomial 
parameterization in Eq. ()91|) and also used a parametrization with an additional x'^ term, 
with free parameter k. In spite of the new parameters, all acceptable results (i.e. with the 
values of comparable to our best fit solution) extracted from fits to data were very similar 
to the ones obtained with parametrization (j91|) . We emphasize that the behaviour of the 
function 6f2{x) for x < 0.70 is well constrained by data and only small variations on both 
the shape and the position of the zero Xq are allowed. This observation in turn results in 
reduced systematic uncertainties of the model. 




FIG. 4: Off-shell function 5f2{x) and effective cross-section axiQ^) corresponding to the parameters 
from Table mi The curves show the size of the uncertainty bands (itlcr ), including both statistical 
and systematic (Table HVj) uncertainties. The effect of different functional forms is also included, 
as explained in Sec. 

For the coherent nuclear effects, we varied (Jq within the uncertainty estimated by av- 
eraging over p^, uj and mesons, ±3 mb f Section IV E|) . As explained in Sec lVGll this 
parameter is strongly correlated with ax and Qq. Similarly, we varied ai within the la 
allowed range fSection lVG We also tried to change the exponent controlling depen- 
dence in Eq.fl^. We obtained almost equally good fits with the dipole and monopole forms 
in Eq.(jH21)- The monopole form of Eq. ()92|) had lower x^/d.o.f. for the overall data set. 

Figure m shows the off-shell function Sf2{x) and the effective cross-section axiQ'^) ob- 
tained in Sec. IV (1^ together with the corresponding total uncertainty bands (including 
both statistical and systematic uncertainties). We comment that in our analysis the off- 
shell correction is treated as the first order correction in the parameter v = {p^ — M^)/M^. 
At large a; > 0.7 the off-shell correction can be as large as 25% for heavy nuclei (see also 
Fig. indicating that higher-order terms in v might not be negligible. This can also be a 
source of systematic uncertainty. However, going beyond the first order in v requires the 
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consideration of higher-order relativistic corrections to nuclear wave and spectral functions 
(see discussion in Sec. lIVvO]) that would also affect the treatment of "standard" FMB effect. 
This goes beyond the scope of the present analysis. We also note that our phenomenological 
function 6/2 is extracted from nuclear data and hence it effectively incorporates additional 
contributions from missing terms. 

The final results are dominated by systematic uncertainties, as can be seen from Tables ITTl 
and II VI However, we note that the magnitude of systematic uncertainties is constrained by 
the value of oi a global fit to data, as explained above. Therefore, the availability of more 
accurate (or with wider coverage of kinematics and nuclei) experimental measurements would 
significantly improve our results. 



VI. DISCUSSION 

We now discuss the results obtained from our fit to nuclear data in Sec. IV Gl In Sec. lVI Al 
we address the problem of the normalization of the nuclear valence quark distribution and 
Sec. IVIBI is focused on the implication of this constraint for our analysis. Section IVIDI is 
devoted to the and A dependence of nuclear effects predicted by our model. In Section 
IVIEI we discuss nuclear effects on the deuteron structure functions. 



A. Nuclear valence quark number 

It is instructive to study the contributions due to different nuclear effects to the normal- 
ization of valence quark distribution in a nucleus. Common wisdom is that this quantity 
should not be corrected by nuclear effects since it counts the baryon number of the system. 
Therefore, it is important to verify if different nuclear effects cancel out in the normaliza- 
tion. In the impulse approximation, i.e. if no shadowing and off-shell effects are taken into 
account, the cancellation of nuclear binding and Fermi motion effects in the normalization 
is explicit and it is guaranteed by the normalization of nucleon distribution function ()3U|) 
to the number of nucleons. It should be also noted that nuclear pions do not contribute to 
nuclear valence distribution. In the presence of off-shell (OS) and nuclear shadowing (NS) 
effects different contributions to the valence quark normalization per one nucleon can be 
written as 

iVvai/A = / q..i/A{x) = iVvai/JV + ^N^f + 6Nll (95) 
Jo 

where N^ai/N = 3 is the number of valence quarks in the nucleon and 

SNZ^ = {v) f dx {qr/N{x)5Ux) - qi/N{x)5f,{x)) , (96) 

6N^l= [ dxg,al/7v(x)57^val(x). (97) 
Jo 

Here gvai/Af = u — u + d — d is the nucleon valence quark distribution, 6fq and 6fg are off- 
shell correction functions for quark and antiquark distributions, 67l^i^i{x) is the shadowing 
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correction to the valence quark distribution and {v) = (jP'—M'^) /M^ is the bound nucleon 
off-shellness averaged over nuclear spectral function (for more details see Sec. lVIlX]) .^" 

In general, the off-shell corrections 5fq and 5fq could be different. Since we phenomeno- 
logically extract the off-shell correction 5/2 from a study of IZ2 data it is difficult to unam- 
biguously disentangle off-shell effects for quark and antiquark distributions. The analysis of 
additional data from either Drell-Yan production or neutrino scattering would be therefore 
important. While we defer a detailed analysis of the existing Drell-Yan data from nuclear 
targets 3| to a future publication, no sensitive neutrino data about nuclear effects on struc- 
ture functions are currently available (see also the discussion in Sec lVIIBj) . In this paper 
we rather try to use simple considerations on the nuclear valence quark number in order 
to test the hypothesis of a single universal off-shell correction for all partons in the bound 
nucleon against the case of different corrections 5fq and 5fq. For this purpose it is enough 
to focus on the high region, where we can use Eq.(|T7j). 




FIG. 5: Relative off-shell {^N^^/Z) and nuclear shadowing {5N^^/?)) corrections to the normaliza- 
tion of the valence quark distribution for ^^Fe and ^*^^Pb nuclei (left and right panel, respectively) 
computed as described in Sec. I VI Al The solid curve shows the sum of the off-shell and the nuclear 
shadowing corrections. 



Let us first assume Sfq{x) = 6fq{x). From Eq.(jl7j) we then also have 5/2(2;) = 6fg{x), that 
implies that we have a universal off-shell function for both quark (valence) and antiquark 
(sea) distributions. We evaluate SN^^ by Eq. as a function ofQ^ using the parameters 
of Sf{x) from Table ITTl and the nucleon valence distribution of [39|. The results for iron 
and lead are reported in Fig. indicating a positive off-shell correction of about 1.5 — 2% 
that decreases with Q'^. We then compute the shadowing correction SN^^ by Eq. ()106b|) 
using the effective cross section extracted from our fits (see Sec. IV E() . The results are shown 
in Fig. El It is important to observe that SN^^ is negative and there is large cancellation 
between off-shell and shadowing effects in the normalization over a wide range of Q^. 

Let us now test a different hypothesis, namely no off-shell effect in the sea of bound 
nucleon 6fq = 0. From Eq.(|17j) we then have for the isoscalar nucleon {u + d)6fq = ^-^2(5/2. 



Note that the normahzation of valence quark distribution is not affected by the order of perturbation 
theory analysis and therefore Ea. H95|) holds to any order in as, unlike the Gross-Llewellyn-Smith sum 
rule |9J| which is corrected by both perturbative |95| and nonperturbative effects. 
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In this case the off-shell correction to valence quark number is dominated by the small x 
region and even becomes divergent. Therefore, this assumption leads to unphysical results 
and we have to rule out this case. 

B. Normalization constraints 

In the following we will favor the assumption of a single universal off-shell function 6f{x), 
according to the discussion of the previous Section. This is supported by the existing 7^2 
data we used to extract the phenomenological off-shell function. The universality of Sf{x) 
should be further verified with both Drell-Yan data and future precise neutrino data. In our 
analysis we use the normalization condition in order to fix parameters of the function Sf{x), 
in particular the parameter Xi. As explained in Sec. IV within all possible values of Xi 
providing comparable descriptions of data (x^) we selected the one minimizing the overall 
correction 5A^°f + SN^^f. 

The function 6f{x) measures the change in the quark-gluon structure of the nucleon in 
nuclear environment. This function is not accessible in experiments with isolated proton 
and/or neutron but can generally be probed in nuclear reactions. The results described in 
Sec. IV G 21 demonstrate that inclusive DIS data have a good sensitivity to off-shell effects, 
allowing a precise determination of this correction. 

The phenomenological cross section in Eq. (j92j) effectively incorporates contributions to 
structure functions due to all twists since it is extracted from data. Higher twists are known 
to be important at low and intermediate and for this reason we should not expect an 
exact cancellation between SN^^ and SN^^ calculated with phenomenological cross section. 
Nevertheless, we observe from Fig. El that the cancellation becomes more accurate at higher 

indicating transition to the leading twist regime. In particular, the exact cancellation 
takes place at ~ 15 GeV^. We performed similar calculation for several nuclei and we 
thus verified that this effect is independent of the choice of the nucleus. 

It should be noted that the nuclear data available in the shadowing region are limited 
to relatively low and for this reason the phenomenological cross section (|92p is not 
constrained at high Q^. In this work we evaluate the effective cross section at high by 
treating the condition SN^^ + SN^^ = as an equation on the cross section. We solve this 
equation numerically using the off-shell function 5/2 (x) from Sec. lVG'2l The resulting cross 
section is presented in Fig. [Altogether with phenomenological cross section extracted from 
our fits. In this paper we use the following simple model for the effective cross section. For 

below the crossing point in Fig. IHlwe use phenomenological cross section (jU^ extracted 
form the fits, and for higher we use the cross section calculated from the normalization 
condition. The difference between the two curves in Fig. El below the crossing point is 
attributed to high-twist effects. 

The function 6f{x) is positive for x < xi (see Fig. H}. This implies a negative off- 
shell correction to the structure functions at small Bjorken x because the offshellness u of a 
bound nucleon is negative. Thus the off-shell correction at small x appears as a leading twist 
shadowing correction. Therefore, in this region there is a certain interplay between nuclear 
effects due to coherent nuclear interactions and off-shell effect. In the region xi < x < xq 



^ Wc obtain 5N°^J3 w —0.5 for iron if wc cut off the contribution of the region a; < 10 ^. Changing the 
lower hmit to x ~ 10^^ increases the magnitude of this correction by about factor of 2. 
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the function Sf{x) is negative that provides an enhancement of bound nucleon structure 
functions. Thus in our approach the antishadowing at x ~ 0.1 is hnked to off-shell effects. 
It is important to note that for the valence distributions there is additional antishadowing 
mechanism due to coherent nuclear interactions. Indeed, the presence of substantial real 
part in the C-odd channel (oa = 1) results in the constructive interference of multiple 
scattering interactions at x ~ 0.1 for valence distributions as will be discussed in more 
detail in Sec. lVllAI 





Phenomenological cross section 

— .„,_^ Effective LT cross section 















1 10 

e^, GeV^ 



FIG. 6: Phenomenological cross section ctt extracted from our fits (solid curve) and effective LT 
cross section (dashed curves) computed for iron and lead nuclei as described in Sec. lVI Al 



C. OfT-shell effect and modification of the nucleon size in nuclei 

From our analysis we obtain a positive off-shell correction at large x > xq. Since the 
virtuality — of the bound nucleon is negative this leads to the suppression of valence 
distribution in the bound nucleon at large x. In order to give a qualitative interpretation 
to this result we consider a simple model of the valence distribution in the nucleon and we 
argue that the behavior of 5/2 at large x observed in data can be related to the increase of 
the nucleon core radius in nuclear environment. 

Let us consider the valence quark distribution in terms of spectral representation Eq. pHj) . 
We will consider a sim ple model in which the spectrum of spectator states is approximated 
by a single mass s [isl ISOj 

/^,/iV = 5(s-s)$(t,p2), ^gg) 

where the function $(t,p^) describes the distribution of valence quarks over t = k"^ in 
the nucleon with the invariant mass p^. For the on-shell nucleon the distribution $(t) is 
characterized by a scale A^. In configuration space this scale should be related to the size of 
the valence quark confinement region ~ (the nucleon core radius). From dimensional 
analysis one can write $(t) = Ct,A~^0(t/A^) where and C„ are dimensionless profile 
function and normalization constant. We found that a simple pole model (j){z) = (1 — z)~"' 
results in a reasonable description of the nucleon valence distribution at large x and high 
Q^. In particular, we obtain a reasonable fit to valence distribution of [s^ at = 15 GeV^ 
and X > 0.2 by taking s = 2.1 GeV^ A^ = 1.2 GeV^ and n = 4.4. 



44 



In order to model off-shell dependence of parton distributions we assume that the nor- 
malization constant C„ and the scale A„ become functions of while the profile function 
and the average mass of spectator states s do not change off-shell. We use Eq.ipHj) in 
order to calculate the off-shell modification of the quark distribution 5fq. The (5-function in 
Eq.(pS|) allows us to integrate over the spectrum of residual system. Inspecting the result- 
ing expression we observe a relation between the derivatives of the quark distribution with 
respect to x and p^. After some algebra we obtain 

^fq = c+ [lngvai(a;)]'a;(l - x)h{x), (99) 

where c = dlnC^/ dlnp^ and A = d\nK^/ dlnp^ taken at = M^. It should be noted that 
Eq.()99p is independent of the specific choice of the profile function 0. 

We use Eq.(jnni) in order to reproduce phenomenological function 6/2 at large x. In 
particular, we fix the parameters c and A in order to reproduce the zero of 5/2 at large x 
(xq) and the slope Sf2{xo). Using s = 2.1 GeV^ we obtain A = 1.03 and c = —2.31. The 
function Sfq{x) by Eq.(|99|) is shown in FigjTI together with the phenomenological function 
5/2 (x). One observes that this simple model agrees with phenomenology at large x but not 
at small x at which effect of the nucleon sea is important. 

The positive sign of the parameter A suggests that the scale parameter decreases in 
nuclear environment since p"^ < for bound nucleon. This in turn indicates the increase 
in the nucleon core Vc in nuclear environment ("swelling" of bound nucleon). In order to 
quantitatively estimate this effect we consider the relative change in the nucleon radius 
StcItc. We have drdvc ~ —\5K^/Is^. The relative change in the scale A^, can be estimated 
as (5A^/A^ = \{p^ — M'^) /JvP, where averaging is taken over bound nucleons. We evaluate 
this quantity using our model spectral function for iron and obtain ~ 9% increase in Tc- 

To conclude this Section we remark that the swelling of bound nucleons was discussed 
in the context of quenching of nuclear longitudinal response function in [t^. The change 
of confinment scale in nuclei in terms of a different approach was discussed in the context 
of the EMC effect in jl^, [t^I • The swelling effect was experimentally constrained to < 30% 
from the analysis of Coulomb sum in [t^. 



D. The and A dependence 

In order to illustrate the dependence of 7^.2 calculated in our approach, in Fig. |S1 we 
plot the ratio 7^2(Fe/7V), where N is isoscalar nucleon {p + n)/2, as a function of x for a 
few fixed Q^. We observe from Fig. |H1 significant variations of the ratio 7^2 with at small 
X < 0.1 and large x > 0.65. This dependence can be attributed to several effects. In 
the nuclear shadowing region at small x the dependence of the ratio 7^2 is due to the 
corresponding dependence of effective cross section ax (see Eq.l^^ and Fig. H)). It must be 
also noted that in the region of x between 0.01 and 0.1 the dependence of 7^.2 is affected 
by the dependence of longitudinal correlation length l/ki (see Section ITV Bl and Eqs. fl^ 
and (IS21))- For 0.1 < x < 0.65 the dependence of 7^2 is negligible. At large x the 
dependence is due to the target mass correction effect by Eq. ()11|) in convolution equations. 

In Fig. El we compare the NMC data on dependence of the ratio 7^2(Sn/C) with our 
calculations. We observe an overall good agreement between data and model calculations 



45 



-1.5 I ' ' ' ' ' ' ^ 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

X 



FIG. 7: Phenomenological off-shell function 5f2{x) (solid) in comparison with 5fq{x) (dashed) 
computed using Eq.ljHHI) as described in text. 
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FIG. 8: Our predictions for the ratio of ^^Fe and isoscalar nucleon structure functions calculated 
for = 5, 10, 20 GeV^. The calculation takes into account the non-isoscalarity correction for iron 
by Eq.dini). 



for all values of x within available region of Q^. However, it should be remarked that 
available data on dependence of nuclear effects are still too scarce to make thorough 
phenomenological studies of this effect. In particular, the correlation between x and for 
fixed target experiments and the lack of information about the distributions of data in 
each of the x bins used (typically only the average is provided) can potentially bias the 
calculations where a significant dependence is expected. 

Figures ^ and |21 show that the model reproduces correctly the ratios 7^2 over a wide 
range of nuclei and kinematic regions. The A dependence of the ratio 7^2 is illustrated 
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FIG. 9: The dependence of the ratio 7^2 (Sn/C) for different values of x as measured by the NMC 
The curves with open squares show the corresponding model calculations. For data points the 
error bars correspond to the sum in quadrature of statistical and systematic uncertainties, while 
the normalization uncertainty is not shown. 



in Fig. ^1 for a few fixed values of x. At small x the A dependence is related to the 
multiple scattering coefficients in Eqs. ()62j) and ()67p. through the nucleon number density 
distributions. The increase in the nuclear shadowing effect with A has "geometrical" origin 
and can be attributed to the rising size of heavy nuclei. At large x the A dependence of the 
ratio 7^2 is determined by the corresponding dependence of parameters of nuclear spectral 
function. The slope of 7^2 as a function of x at intermediate x = 0.5-0.6 increases with A 
because of the corresponding increase in the average separation and kinetic energy of bound 
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nucleons. 



It is interesting to note that at x ~ 0.3 the 



ratio 7^2 depends on neither A nor Q^. 
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FIG. 10: The ratio 7^2(A/D) as a function of A for different values of x. The open squares show 
the model calculations corresponding to the average of the data points. For data points the 
error bars correspond to the sum in quadrature of statistical and systematic uncertainties, while 
the normalization uncertainty is not shown. Two curves calculated at constant values are also 
shown for comparison. 



E. Nuclear effects in Deuterium 



Understanding of nuclear effects in deuterium is an important issue since deuterium 
data are often used as the source of information on the neutron structure functions. As 
explained in Sec. IV Gi the determination of d and u parton distributions is sensitive to 
nuclear corrections to deuterium data f Section IIII B|) . In this Section we apply our model 
with the parameters fixed from fit to data from heavy nuclei (see Table IH} in order to 
calculate nuclear modifications in deuterium and compare our predictions with data. We 
take into account nuclear binding, Fermi motion, off-shell, nuclear pion and shadowing 
corrections as explained in Sec. lIVI It should be emphasized that our approach does not 
require any extrapolation from heavy nuclei to deuterium. 

The ratio of the deuteron and the proton structure functions 7^2(D/p) = F2 I ^2 w^is 
measured by the E665 and NMC collaborations 0, 0| in a wide kinematical region of x and 
Q'^. A comparison with these data provide a good test of the applicability of our model to D, 
since these data were not used in the fits described in Sec. lVGl In Fig. ^2 we show the E665 
and NMC data together with the results of our calculations. A good agreement is found 
between data and the model described in this paper. In particular, our prediction of a small 
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FIG. 11: Comparison of E665 and NMC data to our calculations (curve with open squares) for the 
ratio TZ2{x^D/p). For data points the error bars correspond to the sum in quadrature of statistical 
and systematic uncertainties, while the normalization uncertainty is not shown. 



shadowing effect in D seems to be supported by the measured values of 7^.2 (D/p) at small 
values of x. We note that the ratio 7^2 (D/p) also provides a test of the parton distributions 
used in our calculation and in particular of the difference between d and u quark contents. 
This was not the case for all the remaining data listed in Table |l] which were corrected by 
experiments for the neutron excess, thus providing an effective cancellation of PDFs in the 
ratios n2{A'/A) (Section IVUI)- 
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FIG. 12: Comparison of the extrapolations of E139 nuclear data within the nuclear density model 
of Ref.^3| to our calculations (curve with open squares) for the ratio 7^.2(2;, D/N). The error 
bars of the E139 data points correspond to the sum in quadrature of statistical and systematic 
uncertainties. For each x value, the model calculation was performed at the average of the 
experimental points quoted in Ref . [in| . 



49 




FIG. 13: Our predictions for the ratio TZ2{x,D/N) of deuterium to isoscalar nucleon at a fixed 
Q2 = 10 GeV^. The ±la band is also given, including both statistical and systematic uncertainties. 



Unlike the ratio 7^.2 (D/p) the ratio 7^2 (D/-^) = F,p / F2 " cannot be measured directly 
because a free neutron target is not available. The extraction of 7^2(D/^) from SLAC data 
was discussed in RefJlOf in terms of a phenomenological model of the EMC effect in the 
deuterium. In Ref. [lO| the ratio TZ2(D/N) was extracted by extrapolating the measured 
ratios 7^.2 (v4/D) using the nuclear density model of Ref.ji^]. The key assumption was made 
that the quantity TZ2{A/N) — 1 scales as nuclear number density and it was also assumed that 
this ratio is independent of Q^. The values of 7l2(D/N) were given in for x corresponding 
to the X bins of SLAC data. The results are shown in Fig. fT2l together with our calculation 
of the ratio 7l2(D/N) for the same kinematics of the points presented in ^^l-^^ In Fig- ^1 
we show our prediction for the ratio 7^2(D/A^) at fixed = lOGeV^ and the corresponding 
uncertainty band (±lcr), including model systematics. 



VII. APPLICATIONS 



In this Section we apply our results to evaluate nuclear parton distributions (Sec- 
tion IVll A|) and make the predictions of nuclear effects for neutrino structure functions 
(Section EnH). 



The theoretical uncertainties of such extrapolation were not estimated in Q. See also discussion of these 
points in [5^ . 
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A. Nuclear parton distributions 



The parton distributions are process-independent characteristics of the target in high- 
energy processes. Different phenomenologica l ap proaches to the extraction of nuclear PDFs 
(nPDFs) can be found in Refs. [t^ |8^ [sJ, 8^- It should be remarked at this point that 



physics observables are the cross sections and the structure functions, which include con- 
tributions from all twists. The higher- twist terms are generally process-dependent, can be 
essential in the region of relatively low Q^, and, furthermore, can be substantially affected by 
nuclear environment. Therefore, the applicability of the leading twist approximation must 
be considered in comparison with data as well as in any attempt to extract nPDFs. In our 
approach we derive nPDFs from the analysis of nuclear structure functions f Sections IVG 21 
and IVI B|) that allows us to determine both nPDFs and their uncertainties from existing 
data. However, this paper is not aimed at the full nPDF analysis, which will be published 
elsewhere js^]. We rather want to discuss a few different effects which cause modifications 
of nuclear quark distributions. The numerical results shown in this Section were obtained 
using the NNLO proton and neutron PDFs described in Sec. 1111 iJl 



1. Nuclear convolution 

As discussed in Sec. lIVl in the region of high and large x the nuclear structure func- 
tions can be approximated by incoherent contributions from different nuclear constituents 
which can be presented in a convolution form (see Eqs. (j29|) and (j37|) ). The convolution 
formulas look similar for all type of nuclear structure functions suggesting that the convolu- 
tion equations hold for the parton distributions. We denote Qa/Ti.x.Q'^) the distribution of 
quarks of type a in a target T. Then the quark distribution in a nucleus can be written as 

c=p,n,7r 

where the function fc/A{y, v) can be interpreted as the distribution of particles of type c in a 
nucleus over light-cone momentum y and invariant mass (virtuality) v (for bound nucleons 
and nuclear pions see Eqs. (jHU)) and (jHS|) . respectively). The operation f ® q denotes the 
convolution 

f®q= ! '^f[y^v)q^x/y,Q\v). (102) 

J x<y y 

Equations similar to ()101|) can be written for antiquark and gluon distributions in nuclei. 
Note also that the distribution functions are independent of and, therefore, the 
evolution of nuclear PDFs is goverened by the evolution of PDFs of nuclear constituents. 

In view of applications to complex nuclei with different number of protons and neutrons, 
it is usefuU to sort out the contributions to the convolution equation according to isospin. 
Let us consider the isoscalar and isovector quark distributions, qo = u + d and qi = u — 
d. We first address the contributions from bound protons and neutrons to nuclear quark 
distributions. Assuming exact isospin invariance of PDFs in the proton and neutron we have 
simple relations between the isoscalar and the isovector distributions in the proton and the 
neutron 

qo/p{.x, Q^) = go/n(a;, Q^), (103a) 
gi/p(x,Q^) = -gi/„(x,(5^). (103b) 
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Using these relations we observe that the quark distributions with different isospin decouple 
in the convolution equation. In particular, for the isoscalar (go/A) and the isovector (gi/^) 
nuclear quark distributions we have 

qo/A{x, Q^) = Afo® qo/p, (104a) 
qi/A{x, = (Z - iV) /i ® gi/p, (104b) 

where /o and /i are the isoscalar and the isovector nucleon distributions in a nucleus. These 
distributions are given by Eq. (j^^ with the spectral functions Vq and Vi defined by Eq.fH^. 
Note that the distributions /o and fi are normalized to unity. 

Let us now discuss the pion contribution to Eq. ()10H) . Similar to the nucleon case, we 
assume the isospin relations for quark distributions in the pion: go/7r+ = Qo/w- = %/nO and 
qi/T,+ = —qi/TT- and gi/^o = 0. Using these relations we have for the pion correction to the 
isoscalar and isovector nuclear quark distributions 

?om(^> Q^) = Itt/a ® go/TT, (105a) 
Qi/Ai^^ Q^) = U'tt+za - fn-fA) ® qi/n- (105b) 

Here in the first equation fj^/A is the sum of the distributions over all pion states. It should 
be emphasized that in Eqs. ()105|) the pion distributions refer to nuclear pion excess, since 
scattering off virtual pions emitted and absorbed by same nucleon (nucleon pion cloud) are 
accounted in the proton and neutron PDFs. For the calculation of nuclear pion distributions 
in our model see Sec. IV CI 
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FIG. 14: Nuclear effects for the isoscalar and the isovector quark distributions in ^^Fe. The ratios 
TZo and TZi (see text) were calculated for the valence quark distributions at = 20GeV^. Nuclear 
shadowing and pion corrections are also included for the isoscalar distribution. 

The isovector component should vanish in isoscalar nuclei with Z = N}^ However, 
for a generic nucleus with different number of protons and neutrons both the isoscalar 



It should be remarked that this statement apphes to nuclear states with the total nuclear isospin 0. If 
highcr-isospin states are present for a, Z = N nucleus, then the isovector distribution qi/A may be non-zero. 
The discussion of these issues is postponed for future studies. 
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and the isovector distributions are present. Heavy nuclei typically have a small excess of 
neutrons over the protons and the distributions /o and /i are quite different in such nuclei, 
as discussed in Section IIV A 51 For this reason nuclear effects in PDFs depend on isospin. 
In order to illustrate this statement we calculate the ratios TZq = Qn/A (x) I [A Qn/ p(x) ) and 
T^i = 'li/Aix)/[{Z — N)qi/p{x)] for the iron nucleus using the proton PDFs of Ref.j39[. The 
results are shown in Fig.[T31 We note that the full nuclear correction is shown in case of Qq/a, 
i.e. the calculation includes effect of nuclear spectral function, off-shell correction, nuclear 
pion and shadowing effects. However, for the isovector quark distribution qi/A we neglect 
possible nuclear pion and shadowing effects. 



2. Nuclear shadowing 

In this Section we discuss coherent nuclear effects in the context of parton distribu- 
tions. To this end we want to apply the approach discusses in Section IIV B 21 The multiple 
scattering effects are generally different for different PDFs. We specify this statement by 
considering nuclear effects for quark distributions of different C parity, g^^^(x) = g(x)±g(x). 
In order to simplify discussion we consider the isoscalar (anti)quark distributions, q = u + d 
and q = u + d. The C-odd distribution is in fact the valence quark distribution in the target 
q^~'^ = g^ai- The C-even distribution at small x describes the target quark sea. 

In order to bridge between Sec. IIV Bl and the present discussion we recall that the struc- 
ture function Fi in the LT approximation is given by C-even distribution q^~^\ The structure 
function Fi is transverse helicity structure function (to be more precize, the average over 
left- and right-polarized transverse helicity structure function, see Eq.Q). At small x, as 
discussed in Sec. IIV Bl nuclear effects are described by the propagation of hadronic compo- 
nent of virtual boson with the proper helicity state in nuclear environment. Equation ()63|) 
applies in the case of q^~^\ 

Similarly, the structure function F^ in the LT approximation is given by C-odd (valence) 
distribution q^~\ In terms of helicity structure functions this is the asymmetry between left- 
and right-polarized states. Therefore, nuclear corrections to g'-^-* at small x can be described 
by Eq. (j66|) . We have for coherent nuclear corrections to q^^^ and g'-^-* quark distributions 

Sn^+) = ^1a\^) ^ Ke(a^C2^)/ Ima^, (106a) 
^Iv (^) 

^^(_) ^ Sq^ {x) ^ ^2Re{AaaTC^) - Im(Aa a^Cg^)]/ Im Aa, (106b) 
Qn i^) 

where and are given by Eqs. (j62|) and (j67p with effective transverse scattering am- 
plitude ar = {i + aT)o-T/2. Equation ()106b|l determines the nuclear shadowing effect for 
valence quark distribution (57^vai = STl^"\ The amplitude Aa describes the left-right asym- 
metry in the transverse amplitude. In other terms Aa can be interpreted as the difference 
between qN and qN scattering amplitudes [H^. As discussed in Sec. lIVB'2l the correction 
STl^~^ does not depend on the specific value of the cross section asymmetry Aa but does 
depend on a a = Re Aa/ Im Aa. The rate of nuclear effects for both C-even and C-odd dis- 
tributions is determined by transverse amplitude ay. Nuclear shadowing effect for antiquark 
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distributions can readily be derived from Eqs. ()lU6a|) and ()106b|) and we have 



2qN{x) 



(107) 



The results of calculation of nuclear effects for valence quark and antiquark distributions 
are reported in Fig. El The calculations account of the effects of smearing with nuclear 
spectral function (FMB), off-shell corrections (OS), nuclear shadowing (NS), and nuclear 
pion (PI) corrections. The FMB, OS, and PI corrections have been computed as discussed 
in Sec. lVII A II using our model spectral function, pion distribution function and off-shell 
correction described in Sec.|3 The NS correction for valence and sea distributions are 
computed by Eqs. ()106b|) and ()107|) using the parameters of effective scattering amplitude 
derived from our fits (see Sections IV G 21 and IVI B|) . 
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FIG. 15: Nuclear effects for isoscalar valence and sea quark distributions calculated for iron nucleus 
at = 20GeV^ (see text). The left panel displays different contributions to TZvaV- the dot-dashed 
curve if only the effect of nuclear spectral function (FMB) and off-shell (OS) corrections are taken 
into account, the full curve is overall nuclear correction including nuclear shadowing effect (NS). 
The right panel displays similar contributions to T^sca- The full curve also includes the nuclear 
pion effect (PI), which is absent for the valence quark distribution. 

A few remarks are in order. At small x < 0.01 the NS effect for valence quark distri- 
bution is enhanced relative to that for nuclear sea. The underlying reason for that is the 
enhancement of multiple scattering corrections for the cross section asymmetry as discussed 
in Sec. II V B'2l If we keep only the double scattering correction then the ratio 57^vai/(^'^sea 
is given by Eq. (j60|) . The OS correction is negative in this region. However, the combined 
effect of FMB and OS is somewhat different for valence and sea distributions as displayed 
in Fig. El This is attributed to different x dependence of valence and sea in the nucleon 
which affect the result of the averaging with nuclear spectral function. Nevertheless, in 
spite of these differences, the overall nuclear corrections are similar for valence and sea for 
X < 0.01.^^ 



Note that this discussion refers to a high 
nuclear effects change. 



20 GeV . At lower the balance between different 
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One observes that nuclear corrections for valence and sea distributions are different in the 
antishadowing region. The antishadowing effect for valence {i.e. positive nuclear correction) 
is a joint effect of two corrections both of which are positive: (1) the FMB and OS corrections 
and (2) the constructive interference in the multiple scattering effect which is due to a finite 
real part of the effective scattering amplitude in the C-odd channel (for this reason 
shadowing becomes antishadowing, see the left panel of Fig. II 5|). For sea-quark distribution 
in the antishadowing region we observe a cancellation between different effects. In this 
respect we remark that the contribution of the last term in Eq. ()lU7|) becomes increasingly 
important at x > 0.05, because of the ratio qva.\/N{x)/qN{x). This term is negative in this 
region and cancels a positive nuclear pion contribution. As a result the overall nuclear 
correction to antiquark distribution is small for 0.02 < a; < 0.2. Note that this agrees with 
the results of E772 experiment, in which no enhancement of nuclear sea was observed in DY 
nuclear processes [l^ . 

It should be noted that the calculation of the relative nuclear correction for valence 
quark distribution is stable with respect to the choice of the PDF set for entire region of 
X (see also Fig. El for nuclear correction to valence distributions). Nuclear effects for sea 
quarks also depend weakly on the particular choice of PDF for small x. However, at high x 
the calculation of nuclear effects for antiquark distributions has larger uncertainties and the 
result is sensitive to both the shape and the magnitude of the nucleon antiquark distribution 
(note the val/sea ratio in Eq. ()107|) ). 

B. Neutrino interactions with nuclei 

In this Section we calculate nuclear effects for neutrino charged-current structure func- 
tions using the approach developed in the previous sections. The study of neutrino inter- 
actions is particularly interesting to this end since they are flavour sensitive and they are 
strongly influenced by the structure function F3, which is not present in the electromagnetic 
case. We also note that due to the low interaction probability in practice the detection of 
neutrinos always requires heavy nuclear targets. Therefore, the knowledge of nuclear effects 
is crucial for understanding of neutrino cross-sections. 

In order to compute corrections to F2 and related to the averaging with nuclear 
spectral function (FMB and OS effects) we apply Eqs. (jTTj) and ()28|) and use the off-shell 
function 5/2 extracted from the analysis of Sec. IV G^ for both F2 and F3. Nuclear shadow- 
ing/antishadowing corrections are computed as discussed in Sec. lIVBl 

We focus here on the region of relatively high momentum transfer > 5 GeV^ and as- 
sume that coherent nuclear interactions driven by axial current are similar to those of vector 
current at large and that they can be described by the effective amplitude extracted from 
the analysis of Sections IV El to FV G 31 ^^ A detailed study of nuclear effects in ( ant i) neutrino 
interactions including the low region will be the subject of a future publication [s^l- 

We calculate the ratios = F^^^/^AF^^) and = xF^^^/^AxF^^), where N denotes 
the isoscalar nucleon (averaged over proton and neutron), for the most common nuclear 
targets used by recent neutrino experiments: (NOMAD @), ^'^Fe (NuTeV MI- 
NOS 13), ^°Ar (ICARUS 0) and 207pb (OPERA 0, CHORUS 0). Our results for 



Note that the interactions of the axial-vector current at low are essentially different from those of the 
vector current. This region requires a special analysis which goes beyond the scope of the present paper. 
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FIG. 16: Our predictions for ratios of F2 (left plots) and xF^ (right plots) for neutrino scattering on ^^C and ^^Fe and the corresponding 
values on isoscalar nucleon (p+n)/2. The curves are drawn for = 5, 10, 20 GeV^ and take into account the non-isoscalarity correction. 





FIG. 17: Our predictions for ratios of F2 (left plots) and X-F3 (right plots) for neutrino scattering on ^'^Ar and ^'^'^Pb and the corresponding 
values on isoscalar nucleon (p+n)/2. The curves are drawn for = 5, 10, 20 GeV^ and take into account the non-isoscalarity correction. 



both F2 and xF^ are shown in Figures IT^ and IT7I for different values of . We briefly 
comment on the main features that distinguish the nuclear corrections in neutrino DIS from 
the ones in charged-lepton DIS (-F^) . By comparing Fig. |H1 and El we observe that nuclear 
effects for F2 and in the coherent region are similar (note that we restrict the present 
discussion to relatively high Q^). However, at large x nuclear effects for F2 and F2 are 
somewhat different. In particular, we note that > 7^2 region of x ~ 0.6 — 0.8. 

This is because the neutron excess correction is positive for F2, while it is negative for F^. 
From Fig. El and El one can also observe that nuclear effects at large x are similar for neu- 
trino F2 and XF3. However, at small x the nuclear shadowing effect for xF^ is systematically 
larger as follows from Eq. ()60|) . 
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FIG. 18: Comparison of different nuclear effects calculated for neutrino xF^ at = 10 GeV^ for 
^'^''Pb target. The labels on the curves correspond to the effects included in turn: the averaging with 
nuclear spectral function (FMB), off-shell correction (OS), nuclear pion excess (PI) and coherent 
multiple-scattering correction (NS). The calculation takes into account the target mass and the 
neutron excess corrections. Note that xF^, is not corrected for pion excess effect (overlapping 
dashed and dashed-dotted curves). 

Figure El illustrates different nuclear corrections to xF^ for ^°^Pb target computed at 
fixed Q^. The enhancement at intermediate x values is a joined effect of all considered 
nuclear correction (see also Fig. Eland discussion in Sec. lVII A2|) . In the case of F2 the 
"antishadowing" at x ~ 0.1 is due to off-shell and nuclear pion corrections.^^ Note that 
the nuclear pion excess effect can be neglected in the case of xF^, in contrast to the case of 
F2. Indeed, in the isoscalar nucleus the pion correction depends on pion structure functions 
averaged over different pion states and F^ vanishes after such averaging. A small isovector 



See also Fig. O for F^. Note, however, that the neutron excess correction has a different sign for and 
F2 that explains the differences between the magnitude of nuclear effects in Figs. 151 and II 71 
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correction which is proportional to vr"*" — vr" asymmetry in the nuclear pion distribution 
functions (see Eq. fll05|) ) is also neglected. 

The study of xF^ is particularly important since it allows to test the normalization of 
nuclear valence quark number. As discussed in Sections IVl Al and IVIBI the conservation 
of the nuclear valence quark number was used in our analysis in order to test the balance 
between nuclear shadowing and off-shell effects. The valence quark (baryon) number of 
the target is related to the integral of neutrino and antineutrino averaged F3, the Gross- 
Llewellyn-Smith (GLS) sum rule jo^l- We remark, however, that in QCD this relation is not 
exact and only holds in the leading twist and the leading order in as and is corrected by 
both the radiative [o^ and the higher-twist effects. It would be interesting to experimentally 
address the question of nuclear modification of the GLS sum rule. New measurements of 
xF^ from neutrino and antineutrino scattering off different nuclei would help to clarify this 
issue, provided they can reach a precision comparable to the size of the effects we observe 
in our analysis (typically 1%, see Fig. EI). 

We conclude this Section by remarking that in spite of the major interest of neutrinos as 
a probe for nuclear effects, virtually no precise experimental information is available so far in 
DIS region. The only direct measurements of nuclear effects on neutrino DIS cross-sections 
were performed by BEBC |90| (^ONe/D) and CDHSW lU (^^Fe/p). However, these results 
are affected by large statistical and systematic uncertainties. It should be emphasized that 
neutrino DIS provides information complementary to that of the charged-lepton scatter- 
ing and, therefore, the completion of new high-statistics measurements would have a large 
impact on our understanding of nuclear effects. The NOMAD experiment js^ collected 
large neutrino samples on ^^C, ^^Al and ^^Fe targets allowing a study of nuclear effects from 
27A1/^2Q and ^^Ye/^'^C ratios O- The recent NuTeV cross-section data |9J also provide 
information on nuclear effects in ^^Fe. In addition, the CHORUS experiment j93| is extract- 
ing neutrino cross-sections from the interactions collected on ^"^Pb. Table |V] summarizes 
the various (anti) neutrino data samples. 



Experiment 


Targets 


CC statistics (xlO^) 


Data taking 


Reference 


NOMAD 


12c 


L3(0.06) u{u) 


1995-1998 


[85j 






L5(0.07) u{u) 










12.5(0.6) u{u) 






CHORUS 


207pb 


1.3(0.3) 


1998 


[93] 


NuTeV 


56pg 


2.7(1.2) 


1996-1997 


m 



TABLE V: List of recent nuclear data which can be used to study nuclear effects on neutrino 
structure functions. 



VIII. SUMMARY 

We presented a detailed phenomenological study of unpolarized nuclear structure func- 
tions for a wide kinematical region of x and Q^. A general approach was developed which, 
on one side, includes the main nuclear corrections and, on the other side, provides a good 
description of data on nuclear structure functions. We take into account the QCD treatment 
of the nucleon structure functions and address a number of nuclear effects including nuclear 
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shadowing, Fermi motion and nuclear binding, nuclear pions and off-shell corrections to 
bound nucleon structure functions. 

Starting from a relativistic approach in the description of nuclear DIS we then exploited 
the fact that characteristic energy and momentum of bound nucleon are small compared 
to the nucleon mass. This allowed us to compute nuclear corrections in terms of nonrel- 
ativistic nuclear spectral function, the quantity which is well constrained by data at low- 
and intermediate-energy regions. Our analysis suggested that data cannot be quantitatively 
explained in impulse approximation by applying "standard" Fermi motion and nuclear bind- 
ing corrections even at large x. This motivated us to address the off-shell effect in bound 
nucleon structure functions. This correction was parametrized in terms of a few parameters 
which were extracted from data, together with their uncertainties. The effective scattering 
amplitude which determines the magnitude of nuclear shadowing effect was also addressed 
phenomenologically. 

It should be emphasized that the phenomenological parameters of our model refer to 
the nucleon structure and for this reason they are common to all nuclei. We verified this 
hypothesis by extracting them from different subsets of nuclei. Overall, we obtained an 
excellent agreement between our calculations and data by using only three independent 
parameters. 

Our results show that inclusive nuclear DIS data have a good sensitivity to off-shell 
effects, allowing a precise determination of this correction. We also note that the study of 
semi-inclusive nuclear DIS in which the kinematics of the active nucleon can be controlled by 
selecting certain final states would provide additional information on the off-shell effect. The 
off-shell effect is related to the modification of the nucleon structure in nuclear environment. 
This relation was discussed in terms of a simple model in which the off-shell effect at large 
X was linked to the modification of the bound nucleon core radius. We found that the off- 
shell correction derived from our analysis favours the increase in the nucleon core radius in 
nuclear environment. 

We studied in detail the and A dependencies of nuclear corrections. One important 
application was the calculation of nuclear effects for deuterium, which is of primary interest 
for the problem of the extraction of the neutron structure functions. We also applied our 
model to study nuclear valence and sea quark distributions, as well as the flavour (isospin) 
dependence of nuclear effects. 

Another important application was the calculation of nuclear structure functions for neu- 
trino scattering. In the present paper we evaluated nuclear corrections for charged-current 
neutrino structure functions for relatively high Q^, which are relevant for the analysis of 
existing DIS neutrino data. More detailed studies of neutrino and antineutrino interactions 
for both charged-current and neutral-current scattering are planned in future publications. 
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APPENDIX A: INTEGRATION IN NUCLEAR CONVOLUTION 



The integration in convolution formulas is constrained by the requirement that the in- 
variant mass of bound nucleon and the virtual photon is high enough for producing physical 
final states. In particular, for the region of invariant masses of final states larger than a 
given mass Mx the required relation is 

W'^>Ml, (Al) 

where = {p + q)^ and p the four- momentum of the bound nucleon. The threshold of 
inelastic channels corresponds to Mx = M + m^r and by setting Mx = M we take into 
account the elastic channel. Here we discuss in detail the constraints on the integration 
region in the convolution formulas due to Eq. ()Al|) . Note that Eq. ()Al|) is equivalent to 

Po + qo>Ex, (A2) 

where Ex = {M^ + (p + 9)^)^^^ and po = M + e. Using this equation we can write the 
integral over the bound nucleon four-momentum in convolution formulas as 

J d^p e{W^ -Ml) = j dV j de. (A3) 

Ex-qo-M 

This equation should be applied together with the nuclear spectral function and other func- 
tions which enter the convolution formulas. The energy integration in Eq. ()A3|) corresponds 
to the integration over the excitation energies of the residual nucleus. 
We first consider the spectral function 

V{e,p) = 2Tx5{e-ep)n{jp), (A4) 
£p = ^^o-pV(2mo). (A5) 

This is the relevant case for the deuterium, for which Eq = Ed and mo = M (see Eq. (j39|) ). and 
also for the model spectral function Vmf with mg = M^-i the mass of the residual nucleus 
and Eq = —E^^"* the nucleon separation energy averaged over mean-field configurations of 
the residual nucleus (see Eq.flTIj) and the discussion thereafter). 

The energy integration in Eq. ()A3|) can easily be performed and inequality ()A2|) then 
becomes 

qo + M + Ep>Ex. (A6) 

This inequality provides the constraints on the momentum space in Eq. ()A3|) . In order to 
solve it explicitly we chose the coordinate system such that the momentum transfer has only 
longitudinal component q = (go7 0_L, — |g|). Then after some algebra ()A6|) can be written as 
(we retain only the terms linear in ^p) 

pV(2m,)-p,-p, <0, (AT) 



where the notations are 

= M 



Qy M 

mo + qo + M' 



(A8a) 



m. = (A8b) 
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and A = - M^, 7 = \q\/qo, and 72 = 1 + M/qq. 

Inequality ()A7|) is most easily solved in terms of longitudinal and transverse coordinates, 
p = {pj_,Pz)- In this case, the solution to ()A7|) can be written as 



< Pz ^ t^z, 

< pI < T\ 



(A9) 



where = ml + 2m*p* is maximum transverse momentum squared of the bound nucleon 
(for the given kinematical conditions) and correspond to those longitudinal momenta at 
which the left side of (IA7D is 0, 



The momentum integral in Eq. ()A3|) in terms of these variables is 



dV = 7r^(T) 



dpz 



(AlO) 



(All) 



The reqirement > gives the constraint on possible x and in inelastic scattering off 
bound nucleon. 

In spherical coordinates we introduce the azimuthal angle 6 between the 2-axis and the 
direction of the momentum pz = pcosO (here p = \p\). The solution to ()A7|) splits into two 
different regions with respect to the sign of p* and the momentum integral in convolution 
formula can be written as 



d^p 



p_l_(cos 0) 

2-K I dcos6' / dpp^, if > 0, 



-1 







1 p+(cos6») 

2n I dcos6' J dpp^, 

p_ (cos 9) 



(A12) 



if-^<P*<0, 



where p±(cos6') are the values oi p at which the left side of ()A7|) is 

p±(cos 9) = m^, cos 9 ± a/ (m* cos 9Y + 2m^p^, 



(A13) 



and = (2\p^\/m^)^^'^ . The first case in ()A12|) applies if x < 1 as can be readily seen from 
Eqs. fjA8| ). while the last case concerns the regions x ~ 1 and x > 1. 



The equation = determines the maximum possible x which can be achieved in DIS from bound 
nucleon. In application to the deuteron this gives x = 3/2 (neglecting terms and ea/M corrections in 
Ea. ljA8() ^. This is different from the kinematical maximum x = Mo/M k, 2, which corresponds to elastic 
scattering from the deuteron as a whole. We comment that the limit x = 3/2 was derived keeping linear 
terms in e/M. However, the events with such large x are due to high-momentum configurations p ^ M 
in the wave function and, therefore, require fully rclativistic description. 
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We now consider generic spectrum in Eq. ()A3|) . The upper limit of energy integration 
is determined by the threshold separation energy Sth- We recall that in our notations the 
separation energy e = Eq — E^~^, where Eq is the ground state energy of the target nucleus 
and E^~^ is the energy of the residual nucleus including the recoil energy. Therefore, for 
the given recoil momentum Eth = £o ~ P^/(2?7io), where Eq 



En ^ is the difference of 



the ground state energies of the target and the residual nucleus and mo = Ma^i is the mass 
of the residual nucleus. The constraints on the momentum space in Eq. ()A3|) directly follow 
from 

Ex-qo- M < £th. (A14) 



This inequality, written in terms of Eq and mo, is equivalent to ()A6|) . Therefore, the discussion 
of ()A6|) can be taken over ()A14|) . In particular, the solutions to ()A14|) in terms of the 
longitudinal and transverse momentum are given by Eqs. ()A9|) and ()A10|) . The integration 
in ()A3|) in terms of these variables can be written as 



T2 



dpi 



Pz 



de, 

Ex-qo- M 



(A15) 



where the limits of integration are similar to those in Eq. ()All|) and given by ()A8IA9IA10|) . 
The integration in spherical coordinates is 



/ 



1 p+(cOSe) £th 

2 



dp = < 



27ry dcos6'y dpp J de, if > 

-1 Ex-qo-M 

I p+(cOSe) £th 



(A16) 



m* 



271 / dcos^ / dpp^ / de, if — < p* < 0, 



c* p_(cos6») Ex-qo-M 

where the notations are similar to those in Eq. ()A12|) . 



APPENDIX B: MULTIPLE SCATTERING COEFFICIENTS FOR UNIFORM 
NUCLEAR DENSITY 

The magnitude of coherent nuclear effects in the C-even and C-odd structure functions 
F2 and F3 is determined by the terms and (see Eqs. (jH^ and (jU7j) ). These quantities 
can be computed analytically for uniform density distribution with a sharp edge (square well 
model), pyi(r) = pq6{Ra — \r\), which is a reasonable approximation for large nuclei jH^. 
The nuclear radius Ra in this model is related to the r.m.s. nuclear radius as R\ = 
and the central nuclear density is po = ^/{^^a) with A the number of nucleons. The 
coefficients and are 

C^ = ApoRA^r{y), (Bla) 
Ci = AipoRA)\riy), (Bib) 
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where (/j^Y functions of dimensionless and complex variable y = 2i(po/ — kL)RA 



Vl^iy) = [6 - 3y' - 2y' + Q{y - 1) exp(y)] /y\ (B2a) 
^f'iy) = 12 [-4 + y' + y^S + (2 - yf exp(y)] /y'. (B2b) 

If the real part of the amplitude / and ki can be neglected (which is a reasonable approxi- 
mation for X <^ 0.1), then y = Ra/^j with If = {pQa)~^ the mean free path of the particle 
in a nucleus. 
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